Cezar Marcel DOCA ÎNCOVOIEREA BARELOR ESEU FALSIFICABIL Editura Universitatii din Pitesti - 2007 ISBN: 978-973-690-683-1 CUVÂNT ÎNAINTE Principiul falsificationismului - discutat formal pentru prima data de Karl Popper în 1919-1920, reformulat apoi tot de catre el în 1960 si dezvoltat ulterior de Imre Lakatos - afirma ca o lucrare, pentru a fi utila sau chiar si numai stiintifica, trebuie sa fie falsificabila, adica sa poata fi verificata si infirmata Nu confirmarea acesteia este importanta ci infirmarea ei, adica falsificarea ipotezelor pe cale experimentala sau prin observatii În spiritul rationalismului critic, ipotezele si teoriile sunt considerate adevarate doar pâna la prima lor infirmare Aceasta carte nu se substituie vreunui tratat, manual sau curs de specialitate: nici de rezistenta materialelor, nici de analiza matematica, nici de programare Prezenta lucrarea este cu deosebire o culegere de formule si algoritmi de calcul cu aplicabilitate imediata în studiul încovoierilor statice si vibratorii ale barelor continue, omogene si drepte Toate informatiile cuprinse în eseul de fata sunt oferite cititorului sub GNU Free Documentation Licence AUTORUL Institutul de Cercetari Nucleare Pitesti CUPRINS INTRODUCERE 7 ÎNCOVOIEREA STATICA 9 ECUATII, ALGORITMI SI SOLUTII 25 EVALUARI NUMERICE 51 VIBRATII DE ÎNVOVOIERE 69 VIBRATII INDUSE DE CURGEREA PARALELA 85 CONCLUZII 101 BIBLIOGRAFIE 103 ANEXE 105 INTRODUCERE Sub titlul de mai sus ar trebui sa fie enumerate cel putin câteva idei initiale privind elaborarea, respectiv lecturarea celorlalte capitole Asadar: * desi nu este de specialitatea presupusa, autorul s-a aflat de mai multe ori în situatia de a rezolva probleme de încovoiere statica si în regim vibratoriu pentru diferite configuratii de bare (tuburi, grinzi etc ) montate în instalatii tehnologice; * paginile ce urmeaza reprezinta, în ultima instanta, un rezumat al conspectului cuprinzând principalele ipoteze, ecuatii, solutii exacte, formule de evaluare si algoritmi de calcul, utilizate în studiile teoretice (si experimentale) în domeniu; * prin prezentul eseu se ofera cititorului avizat o relativ cuprinzatoare colectie de informatii, abordabila macar si pentru faptul ca toate acestea nu mai trebuie sa fie "redescoperite" prin diferite lucrari dispersate ÎNCOVOIEREA STATICA Literatura de specialitate dedicata Rezistentei materialelor este infinit mai bogata decât ar putea autorul sa parcurga si sa indice ca referinte bibliografice Prezentul capitol cuprinde doar o sumara trecere în revista a principalelor rezultate teoretice privitoare la încovoierea statica a barelor continue, omogene si drepte, unele aspecte fiind reluate si în completarile din Anexele lucrarii Ipotezele mecanicii corpurilor deformabile Studiul deformarii corpurilor solide este fundament pe urmatoarele ipoteze de calcul: * Ipoteza corpurilor continue, omogene si izotrope - modelele teoretice se elaboreaza folosind: 1) functii continue si 2) constante de material având aceleasi valori în orice punct al corpului solid * Ipoteza identitatii proprietatilor mecanice ale elementului infinit mic cu cele ale corpului solid întreg - în fapt, nu se iau în considerare fortele intercristaline * Ipoteza elasticitatii perfecte - sub anumite valori ale eforturilor unitare, deformatiile se anuleaza odata cu disparitia sarcinilor care le-au generat * Ipoteza proportionalitatii dintre eforturi si deformatii - corpul elastic satisface legea lui Hooke, deformatiile supunându-se principiului suprapunerii efectelor * Ipoteza micilor deformatii, cunoscuta si sub denumirea de ipoteza mentinerii dimensiunilor initiale - corpul sufera: 1) deformatii absolute foarte mici în raport cu dimensiunile sale geometrice, respectiv 2) deformatii specifice neglijabile în raport cu unitatea Functie de complexitatea modelului teoretic, calculele se clasifica în: de ordinul I atunci când ecuatiile de echilibru se scriu pentru starea nedeformata, de ordinul II atunci când se utilizeaza o schema deformata dar se accepta ca deformatiile sunt mici si de ordinul III atunci când, renuntându-se la ipoteza micilor deformatii, se au în vedere deformatiile mari * Principiul lui Saint-Vénant - sisteme diferite de forte echivalente static produc efecte apreciabil diferite doar în punctele de aplicatie * Ipoteza lui Bernoulli: sectiunile plane si normale pe axa unei bare înainte de deformare ramân plane si normale pe aceasta axa si dupa deformare * Ipoteza starii naturale a corpului sau ipoteza absentei tensiunilor initiale Eforturi si solicitari Corpul solid se poate deforma sub actiunea unei forte (totale) R sau / si a unui moment / cuplul (total) M Marimile R si M se numesc eforturi si pun în evidenta actiunea reciproca dintre doua sectiuni ale corpului solid Putând avea directii oarecare în spatiu, eforturile R si M se descompun în: componente normale la planul sectiunii si componente continute în planul acestei sectiuni, cele doua tipuri de proiectii producând, individual, solicitari simple În cazul concret al unei bare (continue, omogene si drepte) forta totala R are: * o componenta N, normala la planul sectiunii, numita forta axiala (fiind aplicata pe axa barei) - produce fie solicitarea de întindere, fie solicitarea de compresiune; * o componenta T, numita forta taietoare, continuta în planul sectiunii (perpendicular pe axa barei) - produce solicitarea de taiere / forfecare La rândul sau, momentul / cuplul total M se descompune în: * momentul de rasucire Mt, al carui vector este dirijat pe axa, deci perpendicular pe planul sectiunii - produce solicitarea de rasucire / torsiune; * momentul încovoietor Mi, cu vectorul continut în planul sectiunii - produce solicitarea de încovoiere Actiunea simultana a doua sau mai multe eforturi poate da nastere unei solicitari compuse Ecuatiile fibrei medii deformate pentru o bara supusa la solicitarea simpla de încovoiere Sub actiunea (doar a) momentului încovoietor M(x) (s-a renuntat la indicele inferior i), axa unei bare orizontale, paralela cu axa de coordonate x, devine o curba plana continua, denumita fibra medie deformata Daca du, dv si dw reprezinta deplasarile / deformarile pe cele trei axe carteziene x, y si z, atunci, în planul (vertical) xOz, fibra deformata are raza de curbura: unde prin E s-a indicat modulul de elasticitate al (materialului) barei, iar I(x) reprezinta momentul ei de inertie Produsul EI(x) se numeste rigiditate Convenindu-se ca sensul pozitiv al rotirii d? este cel orar (deci negativ trigonometric), atunci: Deformându-se, bara sufera atât deplasari în lungul sau - "translatia" u, cât si perpendiculare pe axa - sageata w Neglijând deplasarea u în raport cu sageata w, tangenta unghiului ? se poate scrie: În cazul micilor deformatii: si: rezultat care, introdus în (*) conduce la concluzia ca fibrei medii deformate i se asociaza ecuatia diferentiala de ordinul II: Pentru barele static determinate, adica atunci când momentul M(x) poate fi determinat, în mod unic, direct din conditiile de echilibru, ecuatia diferentiala (1) se integreaza de doua ori si se obtine functia ce descrie deplasarea / sageata perpendiculara pe axa barei: Constantele de integrare C0 si C1 se pot determina: * fie cunoscând sageata w(x) si rotirea ?(x) în acelasi punct x = x* , sau în doua puncte diferite ale barei: x = x1 si x = x2; * fie cunoscând sagetile w(x1) ? w(x2) în doua puncte oarecare ale barei Formulele de calcul (ecuatiile) ce înglobeaza valorile w(x) si ?(x) corespunzatoare capetelor barei, ale caror abscise se obisnuieste a se desemna prin: x = 0 si x = L, reprezinta conditiile la limita Egalitatea (1) mai este denumita si ecuatia diferentiala simplificata Ecuatia diferentiala exacta are expresia: Totusi, atât timp cât, în orice punct al barei, patratul tangentei unghiului de rotire ?(x) are valori neglijabile în raport cu unitatea - cerinta specifica, de altfel, majoritatii problemelor de rezistenta din practica inginereasca - în locul ecuatiei exacte se poate folosi, cu rezultate de încredere, direct ecuatia simplificata Pe de alta parte, daca se tine cont si de efectele fortei taietoare T(x) asupra deformarii fibrei medii, atunci ecuatia diferentiala de ordinul II (simplificata) devine: ; unde G este modulul de forfecare, µ este coeficientul lui Poisson, iar A'(x) este aria redusa a sectiunii transversale si se calculeaza cu formula: Marimea S din ultima formula reprezinta momentul static fata de axa neutra, corespunzator partii din sectiunea transversala care tinde sa lunece, iar b este o dimensiune geometrica a sectiunii De exemplu, în cazul unei sectiunii dreptunghiulare de arie A, calculele conduc la rezultatul A' = 5A / 6, iar în cazul unei sectiuni transversale circulare de arie A se obtine A' = 9A / 10 Pentru barele static nedeterminate momentul de încovoiere M(x) nu mai poate fi exprimat din ecuatiile de echilibru, de aceasta data trebuind sa fie rezolvata ecuatia diferentiala de ordinul IV a fibrei medii deformate: unde q(x) este sarcina distribuita pe lungimea barei Cunoscând legile I(x) si q(x), ecuatia (2) se integreaza de patru ori si se obtine functia deplasarii / sagetii: Constantele de integrare C0, C1, C2 si C3 se determina si de aceasta data impunând solutiei sa respecte, de exemplu, conditiile de rezemare a barei, acestea din urma fiind exprimate prin ecuatii de caracterizare a deplasarilor (sageti si rotiri) si eforturilor (momente si forte taietoare) Functia are dimensiunea unui moment, iar functia are dimensiunea unei forte Rezemarea la capete Generalizând imaginea comuna a unui corp solid sprijinindu-se (în câmp gravitational) pe un al doilea, prin rezemare s-ar putea desemna starea de echilibru mecanic într-un anumit punct (reazem), caracterizata prin valori cunoscute pentru: sageti, unghiuri de rotatie, momente încovoietoare sau forte taietoare Astfel, capatul liber al unei bare montata în consola este un reazem caracterizat prin moment încovoietor si forta taietoare nule Renuntând la montajul în consola, capatul neîncastrat al barei din exemplul anterior s-ar putea "rezema", eventual, pe o forta externa, concentrata, actionând de jos în sus Levitatia magnetica (la capetele unei bare) reprezinta un caz de "sprijinire" fara contact mecanic; si exemplele ar putea continua Ecuatiile de caracterizare a reazemelor se numesc conditii omogene, respectiv conditii neomogene, dupa cum valorile indicate pentru: deplasare, rotire, moment încovoietor, sau forta de forfecare sunt nule, respectiv diferite de zero Conditiile ce se refera la modul de sprijinire în abscisele x* = 0, respectiv x* = L caracterizeaza reazemele de capat Reazemele ale caror abscise îndeplinesc conditia 0 0 I(x) = I0 q(x) = q0 ecuatia (4) devine: si are solutia generala: unde: Daca q0 = 0, atunci: Bara fara pat elastic, forta axiala de compresie, cu moment de inertie constant si sarcina distribuita variabila Fiind date conditiile: kp = 0 N > 0 I(x) = I0 q(x) ecuatia (4) devine: si are solutia asigurata de algoritmul: unde: În cazul particular: avem: Daca q0 = 0, atunci: Bara fara pat elastic, forta axiala de compresie, cu moment de inertie variabil si sarcina distribuita constanta Fiind date conditiile: kp = 0 N > 0 I(x) q(x) = q0 ecuatia (4) devine: careia autorul nu a reusit sa-i gaseasca o solutie generala sau un algoritm de calcul asociat Totusi, în cazul particular: avem solutia asigurata de algoritmul: unde: Prin J0(z) si Y0(z) s-au desemnat functiile Bessel de primul si al doilea ordin Amintim faptul ca functiile Bessel Jn(z) si Yn(z) sunt solutii ale ecuatiei diferentiale: Daca q0 = 0, atunci: Bara fara pat elastic, forta axiala de compresie, cu moment de inertie variabil si sarcina distribuita variabila Fiind date conditiile: kp = 0 N > 0 I(x) q(x) ecuatia (4) devine: careia autorul nu a putut sa-i gaseasca o solutia generala sau un algoritm de calcul asociat Totusi, în cazul particular: avem solutia asigurata de algoritmul: unde: Particularizând, în continuare, cu: avem: Daca q0 = 0, atunci: Bara fara pat elastic, forta axiala de întindere, cu moment de inertie constant si sarcina distribuita constanta Fiind date conditiile: kp = 0 N 0 N = 0 I(x) = I0 q(x) = q0 ecuatia (4) devine: si are solutia generala: unde: Daca q0 = 0, atunci: Evident, solutiile anterioare pot fi rescrise înlocuind functiile exponentiale exp(?x) si exp(-?x) prin intermediul functiilor hiperbolice: sh(?x) si ch(?x) Bara pe pat elastic, fara forta axiala, cu moment de inertie constant si sarcina distribuita variabila Fiind date conditiile: kp > 0 N = 0 I(x) = I0 q(x) ecuatia (4) devine: si are solutia asigurata de algoritmul: unde: În cazul particular: avem: Daca q0 = 0, atunci: Evident, solutiile anterioare pot fi rescrise înlocuind functiile exponentiale exp(?x) si exp(-?x) prin intermediul functiilor hiperbolice: sh(?x) si ch(?x) Bara pe pat elastic, fara forta axiala, cu moment de inertie variabil si sarcina distribuita constanta Fiind date conditiile: kp > 0 N = 0 I(x) q(x) = q0 ecuatia (4) devine: careia autorul nu a putut sa-i gaseasca o solutia generala sau un algoritm de calcul asociat Totusi, în cazul particular: si avem solutia asigurata de algoritmul: unde pFq este functia hipergeometrica generalizata: iar G este functia Meijer, definita, pentru un umar real oarecare, r: Reamintim ca functia G(z) este definita: si, pentru valori n întregi, poate fi privita ca o generalizare a functiei factorial pentru argumentul z complex: Bara pe pat elastic, fara forta axiala, cu moment de inertie variabil si sarcina distribuita variabila Fiind date conditiile: kp > 0 N = 0 I (x) q(x) ecuatia (4) devine: careia autorul nu a reusit sa-i gaseasca o solutie generala sau un algoritm de calcul asociat A se vedea si cazul analizat anterior Bara pe pat elastic, forta axiala de compresie, cu moment de inertie constant si sarcina distribuita constanta Fiind date conditiile: kp > 0 N > 0 I(x) = I0 q(x) = q0 ecuatia (4) devine: si are solutia generala: Daca q0 = 0, atunci: Bara pe pat elastic, forta axiala de compresie, cu moment de inertie constant si sarcina distribuita variabila Fiind date conditiile: kp > 0 N > 0 I(x) = I0 q(x) ecuatia (4) devine: si are solutia asigurata de algoritmul: unde: si În cazul particular: avem: Daca q0 = 0, atunci: Bara pe pat elastic, forta axiala de compresie, cu moment de inertie variabil si sarcina distribuita constanta Fiind date conditiile: kp > 0 N > 0 I(x) q(x) = q0 ecuatia (4) devine: careia autorul nu a reusit sa-i gaseasca o solutie generala sau un algoritm de calcul asociat Bara pe pat elastic, forta axiala de compresie, cu moment de inertie variabil si sarcina distribuita variabila Fiind date conditiile: kp > 0 N > 0 I(x) q(x) ecuatia (4) devine: careia autorul nu a reusit sa-i gaseasca o solutie generala sau un algoritm de calcul asociat Bara pe pat elastic, forta axiala de întindere, cu moment de inertie constant si sarcina distribuita constanta Fiind date conditiile: kp > 0 N 0 N 0 N 0 N 1 then k:=1; case tr[i] of 0:begin {capat liber} a[j+2*k+1,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=0; end; 1:begin {simplu rigid} a[j+2*k+1,j+1]:=sinh(kb*x[i]); a[j+2*k+1,j+2]:=cosh(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=q/N; end; 2:begin {simplu elastic} a[j+2*k+1,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=power(kb,3)*E*I*cosh(kb*x[i]) +ke[i]*sinh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*E*I*sinh(kb*x[i]) +ke[i]*cosh(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],2)/2/N; end; 3:begin {incastrare rigida} a[j+2*k+1,j+1]:=sinh(kb*x[i]); a[j+2*k+1,j+2]:=cosh(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*x[i]*x[i]/2/N; a[j+2*k+2,j+1]:=kb*cosh(kb*x[i]); a[j+2*k+2,j+2]:=kb*sinh(kb*x[i]); a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*x[i]/N; end; 4:begin {incastrare elastica} a[j+2*k+1,j+1]:=sinh(kb*x[i]); a[j+2*k+1,j+2]:=cosh(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+1]:=power(kb,3)*E*I*sinh(kb*x[i]) -ce[i]*kb*cosh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*E*I*cosh(kb*x[i]) -ce[i]*kb*sinh(kb*x[i]); a[j+2*k+2,j+3]:=-ce[i]; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=q*ce[i]*x[i]/N-E*I*q/N; end; 5:begin {rezemare generala} a[j+2*k+1,j+1]:=power(kb,3)*E*I*sinh(kb*x[i]) -ce[i]*kb*cosh(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,3)*E*I*cosh(kb*x[i]) -ce[i]*kb*sinh(kb*x[i]); a[j+2*k+1,j+3]:=-ce[i]; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q*ce[i]*x[i]/N-E*I*q/N; a[j+2*k+2,j+1]:=power(kb,3)*E*I*cosh(kb*x[i]) +ke[i]*sinh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*E*I*sinh(kb*x[i]) +ke[i]*cosh(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],2)/2/N; end; 6:begin {forta in capat liber} a[j+2*k+1,j+1]:=power(kb,2)*sinh(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,2)*cosh(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-q/N; a[j+2*k+2,j+1]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-fc[i]/E/I; end; 7:begin {forta in resort elastic} a[j+2*k+1,j+1]:=power(kb,2)*sinh(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,2)*cosh(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-q/N; a[j+2*k+2,j+1]:=ke[i]/E/I*sinh(kb*x[i]) +power(kb,3)*cosh(kb*x[i]); a[j+2*k+2,j+2]:=ke[i]/E/I*cosh(kb*x[i]) +power(kb,3)*sinh(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]/E/I; a[j+2*k+2,j+4]:=ke[i]/E/I; b[j+2*k+2]:=fc[i]/E/I -q*ke[i]*power(x[i],2)/2/E/I/N; end; 8:begin {simplu rigid intermediar} a[j+2*k+1,j+1]:=sinh(kb*x[i]); a[j+2*k+1,j+2]:=cosh(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+5]:=sinh(kb*x[i]); a[j+2*k+2,j+6]:=cosh(kb*x[i]); a[j+2*k+2,j+7]:=x[i]; a[j+2*k+2,j+8]:=1; b[j+2*k+2]:=-q*power(x[i],2)/2/N; a[j+2*k+3,j+1]:=kb*cosh(kb*x[i]); a[j+2*k+3,j+2]:=kb*sinh(kb*x[i]); a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-kb*cosh(kb*x[i]); a[j+2*k+3,j+6]:=-kb*sinh(kb*x[i]); a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=-power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+6]:=-power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 9:begin {simplu elastic intermediar} a[j+2*k+1,j+1]:=E*I*power(kb,3)*cosh(kb*x[i]) +ke[i]*sinh(kb*x[i]); a[j+2*k+1,j+2]:=E*I*power(kb,3)*sinh(kb*x[i]) +ke[i]*cosh(kb*x[i]); a[j+2*k+1,j+3]:=ke[i]*x[i]; a[j+2*k+1,j+4]:=ke[i]; b[j+2*k+1]:=-ke[i]*q*power(x[i],2)/2/N; a[j+2*k+2,j+5]:=E*I*power(kb,3)*cosh(kb*x[i]) +ke[i]*sinh(kb*x[i]); a[j+2*k+2,j+6]:=E*I*power(kb,3)*sinh(kb*x[i]) +ke[i]*cosh(kb*x[i]); a[j+2*k+2,j+7]:=ke[i]*x[i]; a[j+2*k+2,j+8]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],2)/2/N; a[j+2*k+3,j+1]:=kb*cosh(kb*x[i]); a[j+2*k+3,j+2]:=kb*sinh(kb*x[i]); a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-kb*cosh(kb*x[i]); a[j+2*k+3,j+6]:=-kb*sinh(kb*x[i]); a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=-power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+6]:=-power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 10:begin {forta concentrata intermediar} a[j+2*k+1,j+1]:=sinh(kb*x[i]); a[j+2*k+1,j+2]:=cosh(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; a[j+2*k+1,j+5]:=-sinh(kb*x[i]); a[j+2*k+1,j+6]:=-cosh(kb*x[i]); a[j+2*k+1,j+7]:=-x[i]; a[j+2*k+1,j+8]:=-1; b[j+2*k+1]:=0; a[j+2*k+2,j+1]:=kb*cosh(kb*x[i]); a[j+2*k+2,j+2]:=kb*sinh(kb*x[i]); a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; a[j+2*k+2,j+5]:=-kb*cosh(kb*x[i]); a[j+2*k+2,j+6]:=-kb*sinh(kb*x[i]); a[j+2*k+2,j+7]:=-1; a[j+2*k+2,j+8]:=0; b[j+2*k+2]:=0; a[j+2*k+3,j+1]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+3,j+2]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+3,j+3]:=0; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-power(kb,3)*sinh(kb*x[i]); a[j+2*k+3,j+6]:=-power(kb,3)*cosh(kb*x[i]); a[j+2*k+3,j+7]:=-0; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+2]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=power(kb,3)*cosh(kb*x[i]); a[j+2*k+4,j+6]:=power(kb,3)*sinh(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=-fc[i]/E/I; j:=4*(i-1); end; end; end; if N=0 then for i:=1 to Nr do begin if i>1 then k:=1; case tr[i] of 0:begin {capat liber} a[j+2*k+1,j+1]:=x[i]; a[j+2*k+1,j+2]:=1; a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-q*power(x[i],2)/2/E/I; a[j+2*k+2,j+1]:=1; a[j+2*k+2,j+2]:=0; a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*x[i]/E/I; end; 1:begin {simplu rigid} a[j+2*k+1,j+1]:=power(x[i],3)/6; a[j+2*k+1,j+2]:=x[i]*x[i]/2; a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],4)/24/E/I; a[j+2*k+2,j+1]:=x[i]; a[j+2*k+2,j+2]:=1; a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*power(x[i],2)/2/E/I; end; 2:begin {simplu elastic} a[j+2*k+1,j+1]:=x[i]; a[j+2*k+1,j+2]:=1; a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-q*power(x[i],2)/2/E/I; a[j+2*k+2,j+1]:=ke[i]*power(x[i],3)/6+E*I; a[j+2*k+2,j+2]:=ke[i]*power(x[i],2)/2; a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=-q*x[i]-ke[i]*q*power(x[i],4)/24/E/I; end; 3:begin {incastrare rigida} a[j+2*k+1,j+1]:=power(x[i],3)/6; a[j+2*k+1,j+2]:=power(x[i],2)/2; a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],4)/24/E/I; a[j+2*k+2,j+1]:=power(x[i],2)/2; a[j+2*k+2,j+2]:=x[i]; a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*power(x[i],3)/6/E/I; end; 4:begin {inacstrare elastica} a[j+2*k+1,j+1]:=power(x[i],3)/6; a[j+2*k+1,j+2]:=power(x[i],2)/2; a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],4)/24/E/I; a[j+2*k+2,j+1]:=E*I*x[i]-ce[i]*power(x[i],2)/2; a[j+2*k+2,j+2]:=E*I-ce[i]*x[i]; a[j+2*k+2,j+3]:=-ce[i]; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=q*ce[i]*power(x[i],3)/6/E/I -q*power(x[i],2)/2; end; 5:begin {rezemare generala} a[j+2*k+1,j+1]:=E*I*x[i]-ce[i]*power(x[i],2)/2; a[j+2*k+1,j+2]:=E*I-ce[i]*x[i]; a[j+2*k+1,j+3]:=-ce[i]; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q*ce[i]*power(x[i],3)/6/E/I -q*power(x[i],2)/2; a[j+2*k+2,j+1]:=ke[i]*power(x[i],3)/6+E*I; a[j+2*k+2,j+2]:=ke[i]*power(x[i],2)/2; a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=q*x[i]-q*ke[i]*power(x[i],4)/24/E/I; end; 6:begin {forta in capat liber} a[j+2*k+1,j+1]:=1; a[j+2*k+1,j+2]:=0; a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-fc[i]/E/I-q*x[i]/E/I; a[j+2*k+2,j+1]:=x[i]; a[j+2*k+2,j+2]:=1; a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*power(x[i],2)/2/E/I; end; 7:begin {forta in resort elastic} a[j+2*k+1,j+1]:=x[i]; a[j+2*k+1,j+2]:=1; a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=-q*power(x[i],2)/2/E/I; a[j+2*k+2,j+1]:=1+ke[i]*power(x[i],3)/6/E/I; a[j+2*k+2,j+2]:=ke[i]*power(x[i],2)/2/E/I; a[j+2*k+2,j+3]:=ke[i]*x[i]/E/I; a[j+2*k+2,j+4]:=ke[i]/E/I; b[j+2*k+2]:=fc[i]/E/I -q*ke[i]*power(x[i],4)/24/E/E/I/I -q*x[i]/E/I; end; 8:begin {simplu rigid intermediar} a[j+2*k+1,j+1]:=power(x[i],3)/6; a[j+2*k+1,j+2]:=power(x[i],2)/2; a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],4)/24/E/I; a[j+2*k+2,j+5]:=power(x[i],3)/6; a[j+2*k+2,j+6]:=power(x[i],2)/2; a[j+2*k+2,j+7]:=x[i]; a[j+2*k+2,j+8]:=1; b[j+2*k+2]:=-q*power(x[i],4)/24/E/I; a[j+2*k+3,j+1]:=power(x[i],2)/2; a[j+2*k+3,j+2]:=x[i]; a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-power(x[i],2)/2; a[j+2*k+3,j+6]:=-x[i]; a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=x[i]; a[j+2*k+4,j+2]:=1; a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=-x[i]; a[j+2*k+4,j+6]:=-1; a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 9:begin {simplu elastic intermediar} a[j+2*k+1,j+1]:=E*I+ke[i]*power(x[i],3)/6; a[j+2*k+1,j+2]:=ke[i]*power(x[i],2)/2; a[j+2*k+1,j+3]:=ke[i]*x[i]; a[j+2*k+1,j+4]:=ke[i]; b[j+2*k+1]:=-ke[i]*q*power(x[i],4)/24/E/I-q*x[i]; a[j+2*k+2,j+5]:=E*I+ke[i]*power(x[i],3)/6; a[j+2*k+2,j+6]:=ke[i]*power(x[i],2)/2; a[j+2*k+2,j+7]:=ke[i]*x[i]; a[j+2*k+2,j+8]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],4)/24/E/I-q*x[i]; a[j+2*k+3,j+1]:=power(x[i],2)/2; a[j+2*k+3,j+2]:=x[i]; a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-power(x[i],2)/2; a[j+2*k+3,j+6]:=-x[i]; a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=x[i]; a[j+2*k+4,j+2]:=1; a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=-x[i]; a[j+2*k+4,j+6]:=-1; a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 10:begin {forta concentrata intermediar} a[j+2*k+1,j+1]:=power(x[i],3)/6; a[j+2*k+1,j+2]:=power(x[i],2)/2; a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; a[j+2*k+1,j+5]:=-power(x[i],3)/6; a[j+2*k+1,j+6]:=-power(x[i],2)/2; a[j+2*k+1,j+7]:=-x[i]; a[j+2*k+1,j+8]:=-1; b[j+2*k+1]:=0; a[j+2*k+2,j+1]:=power(x[i],2)/2; a[j+2*k+2,j+2]:=x[i]; a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; a[j+2*k+2,j+5]:=-power(x[i],2)/2; a[j+2*k+2,j+6]:=-x[i]; a[j+2*k+2,j+7]:=-1; a[j+2*k+2,j+8]:=0; b[j+2*k+2]:=0; a[j+2*k+3,j+1]:=x[i]; a[j+2*k+3,j+2]:=1; a[j+2*k+3,j+3]:=0; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-x[i]; a[j+2*k+3,j+6]:=-1; a[j+2*k+3,j+7]:=0; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=1; a[j+2*k+4,j+2]:=0; a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=-1; a[j+2*k+4,j+6]:=0; a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=-fc[i]/E/I; j:=4*(i-1); end; end; end; if N>0 then for i:=1 to Nr do begin if i>1 then k:=1; case t[i] of 0:begin {capat liber} a[j+2*k+1,j+1]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=0; end; 1:begin {simplu rigid} a[j+2*k+1,j+1]:=sin(kb*x[i]); a[j+2*k+1,j+2]:=cos(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+1]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=q/N; end; 2:begin {simplu elastic} a[j+2*k+1,j+1]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=-power(kb,3)*E*I*cos(kb*x[i]) +ke[i]*sin(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*E*I*sin(kb*x[i]) +ke[i]*cos(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],2)/2/N; end; 3:begin {incastrare rigida} a[j+2*k+1,j+1]:=sin(kb*x[i]); a[j+2*k+1,j+2]:=cos(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+1]:=kb*cos(kb*x[i]); a[j+2*k+2,j+2]:=-kb*sin(kb*x[i]); a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-q*x[i]/N; end; 4:begin {inacstrare elastica} a[j+2*k+1,j+1]:=sin(kb*x[i]); a[j+2*k+1,j+2]:=cos(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+1]:=-power(kb,3)*E*I*sin(kb*x[i]) -ce[i]*kb*cos(kb*x[i]); a[j+2*k+2,j+2]:=-power(kb,3)*E*I*cos(kb*x[i]) +ce[i]*kb*sin(kb*x[i]); a[j+2*k+2,j+3]:=-ce[i]; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=q*ce[i]*x[i]/N-E*I*q/N; end; 5:begin {rezemare generala} a[j+2*k+1,j+1]:=-power(kb,3)*E*I*sin(kb*x[i]) -ce[i]*kb*cos(kb*x[i]); a[j+2*k+1,j+2]:=-power(kb,3)*E*I*cos(kb*x[i]) +ce[i]*kb*sin(kb*x[i]); a[j+2*k+1,j+3]:=-ce[i]; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q*ce[i]*x[i]/N-E*I*q/N; a[j+2*k+2,j+1]:=-power(kb,3)*E*I*cos(kb*x[i]) +ke[i]*sin(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*E*I*sin(kb*x[i]) +ke[i]*cos(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]; a[j+2*k+2,j+4]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*x[i]*x[i]/2/N; end; 6:begin {forta in capat liber} a[j+2*k+1,j+1]:=power(kb,2)*sin(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,2)*cos(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+2,j+2]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+2,j+3]:=0; a[j+2*k+2,j+4]:=0; b[j+2*k+2]:=-fc[i]/E/I; end; 7:begin {forta in resort elastic} a[j+2*k+1,j+1]:=power(kb,2)*sin(kb*x[i]); a[j+2*k+1,j+2]:=power(kb,2)*cos(kb*x[i]); a[j+2*k+1,j+3]:=0; a[j+2*k+1,j+4]:=0; b[j+2*k+1]:=q/N; a[j+2*k+2,j+1]:=ke[i]/E/I*sin(kb*x[i]) -power(kb,3)*cos(kb*x[i]); a[j+2*k+2,j+2]:=ke[i]/E/I*cos(kb*x[i]) +power(kb,3)*sin(kb*x[i]); a[j+2*k+2,j+3]:=ke[i]*x[i]/E/I; a[j+2*k+2,j+4]:=ke[i]/E/I; b[j+2*k+2]:=fc[i]/E/I -q*ke[i]*power(x[i],2)/2/E/I/N; end; 8:begin {simplu rigid intermediar} a[j+2*k+1,j+1]:=sin(kb*x[i]); a[j+2*k+1,j+2]:=cos(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; b[j+2*k+1]:=-q*power(x[i],2)/2/N; a[j+2*k+2,j+5]:=sin(kb*x[i]); a[j+2*k+2,j+6]:=cos(kb*x[i]); a[j+2*k+2,j+7]:=x[i]; a[j+2*k+2,j+8]:=1; b[j+2*k+2]:=-q*power(x[i],2)/2/N; a[j+2*k+3,j+1]:=kb*cos(kb*x[i]); a[j+2*k+3,j+2]:=-kb*sin(kb*x[i]); a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-kb*cos(kb*x[i]); a[j+2*k+3,j+6]:=kb*sin(kb*x[i]); a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=-power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+2]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+6]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 9:begin {simplu elastic intermediar} a[j+2*k+1,j+1]:=-E*I*power(kb,3)*cos(kb*x[i]) +ke[i]*sin(kb*x[i]); a[j+2*k+1,j+2]:=E*I*power(kb,3)*sin(kb*x[i]) +ke[i]*cos(kb*x[i]); a[j+2*k+1,j+3]:=ke[i]*x[i]; a[j+2*k+1,j+4]:=ke[i]; b[j+2*k+1]:=-ke[i]*q*power(x[i],2)/2/N; a[j+2*k+2,j+5]:=-E*I*power(kb,3)*cos(kb*x[i]) +ke[i]*sin(kb*x[i]); a[j+2*k+2,j+6]:=E*I*power(kb,3)*sin(kb*x[i]) +ke[i]*cos(kb*x[i]); a[j+2*k+2,j+7]:=ke[i]*x[i]; a[j+2*k+2,j+8]:=ke[i]; b[j+2*k+2]:=-ke[i]*q*power(x[i],2)/2/N; a[j+2*k+3,j+1]:=kb*cos(kb*x[i]); a[j+2*k+3,j+2]:=-kb*sin(kb*x[i]); a[j+2*k+3,j+3]:=1; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=-kb*cos(kb*x[i]); a[j+2*k+3,j+6]:=kb*sin(kb*x[i]); a[j+2*k+3,j+7]:=-1; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=-power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+2]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+6]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=0; j:=4*(i-1); end; 10:begin {forta concentrata intermediar} a[j+2*k+1,j+1]:=sin(kb*x[i]); a[j+2*k+1,j+2]:=cos(kb*x[i]); a[j+2*k+1,j+3]:=x[i]; a[j+2*k+1,j+4]:=1; a[j+2*k+1,j+5]:=-sin(kb*x[i]); a[j+2*k+1,j+6]:=-cos(kb*x[i]); a[j+2*k+1,j+7]:=-x[i]; a[j+2*k+1,j+8]:=-1; b[j+2*k+1]:=0; a[j+2*k+2,j+1]:=kb*cos(kb*x[i]); a[j+2*k+2,j+2]:=-kb*sin(kb*x[i]); a[j+2*k+2,j+3]:=1; a[j+2*k+2,j+4]:=0; a[j+2*k+2,j+5]:=-kb*cos(kb*x[i]); a[j+2*k+2,j+6]:=kb*sin(kb*x[i]); a[j+2*k+2,j+7]:=-1; a[j+2*k+2,j+8]:=0; b[j+2*k+2]:=0; a[j+2*k+3,j+1]:=-power(kb,3)*sin(kb*x[i]); a[j+2*k+3,j+2]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+3,j+3]:=0; a[j+2*k+3,j+4]:=0; a[j+2*k+3,j+5]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+3,j+6]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+3,j+7]:=0; a[j+2*k+3,j+8]:=0; b[j+2*k+3]:=0; a[j+2*k+4,j+1]:=-power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+2]:=power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+3]:=0; a[j+2*k+4,j+4]:=0; a[j+2*k+4,j+5]:=power(kb,3)*cos(kb*x[i]); a[j+2*k+4,j+6]:=-power(kb,3)*sin(kb*x[i]); a[j+2*k+4,j+7]:=0; a[j+2*k+4,j+8]:=0; b[j+2*k+4]:=-fc[i]/E/I; j:=4*(i-1); end; end; end; end; Evident, aceasta procedura poate fi completata oricând cu instructiuni de atribuire specifice si altor tipuri de rezeme discutate în capitolul precedent Valoarea determinantului asociat matricei ||a|| se poate calcula, de exemplu, folosind metoda lui Gauss, implementata în corpul functiei mai jos listate: function gauss(n:integer):extended; label 1,2; var i,j,k:integer; p,g:extended; begin p:=1; for i:=1 to n-1 do begin if a[i,i]=0 then begin for j:=i+1 to n do begin if a[j,i]<>0 then begin for k:=1 to n do begin g:=a[i,k]; a[i,k]:=a[j,k]; a[j,k]:=g; end; p:=-p; goto 1; end; end; g:=0; goto 2; end; 1: for j:=i+1 to n do begin if a[j,i]<>0 then for k:=1 to n do a[j,k]:=a[j,k]-a[j,i]*a[i,k]/a[i,i] end; end; g:=1; for i:=1 to n do g:=g*a[i,i]; g:=g*p; 2: gauss:=g; end; VIBRATII DE ÎNCOVOIERE Daca eforturile externe ce actioneaza asupra unei bare se modifica în timp, atunci si deformatiile suferite de aceasta vor depinde de o a doua variabila independenta, t, exprimând faptul ca deplasarile du, dv si dw nu ramân constante În aceste conditii, pentru aceeasi ipoteza du = 0, ecuatia (4), de exemplu, se va rescrie sub forma ecuatiei de miscare: ce s-ar mai putea completa cu fortele de frecare: externa, si interna, , ambele presupuse a fi de natura vâscoasa si caracterizate prin coeficientii asociati cv, respectiv µ Sarcina distribuita q(x,t) va cuprinde, de aceasta data, atât fortele externe excitatoare p(x,t) cât si forta de inertie ; prin ? s-a desemnat densitatea barei iar A este aria sectiunii sale transversale În cazul ideal al unei bare de rigiditate EI constanta, fara frecari interne (adica µ = 0) si în lipsa unor forte de frecare externe (cv = 0, în vid), ecuatia de miscare se simplifica si devine: iar functia w(x,t) trebuie sa satisfaca, pe lânga conditiile la limita cu care ne-am obisnuit deja, si conditiile initiale, precizate pentru t = 0: Moduri proprii de vibratii Din motive ce tin exclusiv de pastrarea formalismului consacrat al lucrarilor de specialitate, în cele ce urmeaza ne vom referi la deformari / deplasari / încovoieri dv(x,t) ce au loc pe directia y Conform definitiilor, modurile proprii (sau naturale) de vibratii ale unei bare ideale sunt caracterizate prin solutiile ecuatiei de miscare în vid, obtinute în lipsa fortelor excitatoare, adica: Utilizând metoda separarii variabilelor, pentru ecuatia (5) se cauta o solutie de forma: si avem: rezultat ce presupune rezolvarea simultana a egalitatilor independente: Concret, plecând de la ecuatia: rescrisa: se cauta solutia generala: care conduce la: de unde rezulta Cu ajutorul solutiilor imaginare: ; functia T(t) devine: sau, exprimând exponentialele exp(i?t) si exp(-i?t) prin intermediul functiilor trigonometrice sin(?t) si cos(?t): ceea ce înseamna ca bara studiata executa oscilatii armonice cu pulsatia proprie: În mod similar, plecând de la ecuatia: avem, consecutiv: ; ; ; astfel încât functia proprie Y(x) capata forma: Exprimând, acum, exponentialele exp(i?t) si exp(-i?t) prin intermediul functiilor trigonometrice sin(?t) si cos(?t), respectiv exponentialele exp(?t) si exp(-?t) prin intermediul functiilor hiperbolice sh(?t) si ch(?t), se regaseste expresia consacrata: Functia proprie împreuna cu pulsatia proprie definesc ceea ce se numeste un mod propriu si fac obiectul studiilor de analiza modala Necunoscutele reprezentând constantele de integrare C1, C2, C3 si C4, respectiv parametrul modal a se determina, si de aceasta data, din conditiile la limita, iar T0 si f din conditiile initiale De exemplu, daca bara de lungime L este încastrata rigid la ambele capete si nu are alte reazeme intermediare, atunci din conditiile la limita se obtine sistemul compatibil nedeterminat de patru ecuatii liniare: ce va admite solutii nebanale daca si numai daca determinantul sau se anuleaza, adica: Ultima egalitate reprezinta ecuatia pulsatiilor proprii si, adusa la forma simplificata: are solutiile: ; ; (n > 2) carora le corespund frecventele proprii: ; (n > 2) Se observa ca exista un numar infinit de moduri de vibratie, desemnate prin valorile naturale ale indicelui n Alegând acum valorile C1n drept variabile independente, atunci pentru celelalte necunoscute C2n, C3n si C4n se obtin solutiile: si forma fibrei deformate în modul n de vibratie va fi descrisa de functia proprie: Evident, pentru oscilatiile armonice vom avea: ; si modul n de vibratie va fi caracterizat, în final, prin: cu constantele An si fn determinate din conditiile initiale (t = 0) În mod similar se poate demonstra ca daca aceeasi bara s-ar sprijini la capete pe reazeme simplu rigide, atunci modurile sale proprii de vibratii ar fi caracterizate cu ajutorul relatiilor: ; Comparând cele doua exemple de încovoiere la vibratii se poate remarca existenta unei dependente directe, explicite, atât a formei fibrei deformate cât si a frecventelor proprii de natura reazemelor (de capat) Vibratiile unei bare reale în mediu fluid vâscos Asa cum am subliniat deja, modurile proprii caracterizeaza situatia ideala a vibratiilor libere (fara forta excitatoare) în vid, ale unei bare fara frecari interne În vederea studierii vibratiilor libere într-un mediu fluid oarecare (vâscos) va trebui sa se tina seama, pe lânga forta de frecare externa, si de faptul ca, oscilând, bara va trebui sa deplaseze (sa disloce) o masa aditionala de fluid, ma (pe unitatea de lungime) Daca se iau în calcul si fortele de frecare interna si externa amintite la începutul acestui capitol, ecuatia de miscare capata expresia: a carei rezolvare permite obtinerea formulei de calcul pentru coeficientul total de amortizare, n: respectiv obtinerea formulei de calcul pentru pulsatiile libere în fluid: Prezentei fortelor de frecare determina executarea de catre bara studiata a unei miscari periodice amortizate, caracterizata prin raportul de amortizare: si, în cazul amortizarii subcritice (? 2 < 1), prin pseudopulsatia: ; Cu ajutorul acestor marimi, solutia ecuatiei de miscare se va scrie: unde functia fibrei deformate Y(ax) se determina din conditiile de rezemare, iar constantele de integrare To si f din conditiile initiale (t = 0) Avem, prin urmare, o oscilatie libera amortizata Masa aditionala Prin masa aditionala (sau masa hidrodinamica) se întelege, în general, masa de fluid deplasata de bara atunci când aceasta din urma executa o miscare vibratorie de încovoiere Definitia de mai sus se poate folosi ca atare numai atunci când avem o singura bara imersata într-un volum infinit de fluid În cazul unui numar oarecare de bare (paralele) imersate într-un fluid (apa) cu volum limitat, masa aditionala este afectata atât de miscarea celorlalte bare precum si de pozitia acestora fata de frontiera fluidului Determinari minutioase ale masei aditionale ma pentru o bara cilindrica de masa mc se pot face masurând frecventele de rezonanta în aer faer si în apa fapa, prin prelucrarea carora obtine formula de calcul: Plecând de la aceste rezultate si tinând cont de definitia data anterior, se poate considera ca masa de fluid deplasata este proportionala cu volumul cilindrului, adica : unde ?f este densitatea fluidului, d diametrul cilindrului (barei) si cm un factor, numit adesea si coeficient de masa aditionala Coeficientul de masa aditionala tine cont de conditiile de frontiera În cazul unui singur cilindru de diametru d, aflat în interiorul unui tub de diametru D, în urma rezolvarii ecuatiilor Navier-Stokes pentru miscarea unidimensionala (cei doi cilindri se considera ca având lungimi infinite), se obtine: Rezolvând aceleasi ecuatii de miscare Navier-Stokes dar de data aceasta în cazul miscarii bidimensionale a unui singur cilindru pozitionat concentric in interiorul unui tub încarcat cu fluid, avem: cu: ; unde Jm si Ym sunt functiile Bessel de ordinul I si II, ? = 2p? pulsatia barei cilindrice iar c este viteza sunetului în fluidul respectiv Formula de calcul se complica semnificativ daca cei doi cilindrii sunt excentrici, si înca si mai mult daca se iau în considerare lungimile lor finite Vibratii întretinute Vibratiile libere de încovoiere ale unei bare drepte pot fi declansate, de exemplu, printr-un soc mecanic (exterior) Asa cum s-a vazut în paragrafele anterioare, prezenta fortelor de frecare interna si externa va determin executarea de catre bara studiata a unor vibratii de încovoiere amortizate Aceasta înseamna ca la un moment dat miscarea va înceta Prioritare, însa, din punct de vedere al studiilor tehnologice, ramân situatiile când aceste vibratii de încovoiere sunt întretinute prin actiunea externa a asa-numitelor forte excitatoare Doua tipuri mai des întâlnite de forte excitatoare ar fi cele armonice, respectiv cele aleatoare Forta excitatoare armonica Vom presupune ca avem de-a face cu forta excitatoare sinusoidala: ce ne permite sa rescriem ecuatia de miscare sub forma: Solutia ecuatiei de miscare fara forta excitatoare reprezinta, asa cum s-a prezentat anterior, o miscare libera amortizata descrisa prin formula: Conform teoriei generale a sistemelor oscilante cu amortizare vâscoasa, pentru ecuatia (6) se cauta o solutie particulara de forma: astfel încât solutia generala a ecuatiei de miscare (6) va fi: Se observa ca, datorita amortizarii, vibratiile proprii se anuleaza exponential astfel încât, dupa trecerea fazei tranzitorii, se poate lua în considerare numai solutia stationara corespunzatoare vibratiilor întretinute Cu alte cuvinte, sub actiunea unei forte excitatoare periodice, bara va executa vibratii de încovoiere având pulsatia egala cu pulsatia fortei respective În acelasi timp functia Yparticular(x) trebuie sa satisfaca atât ecuatia de miscare (6) cât si conditiile la limita în reazeme Acestea vor fi implicit îndeplinite daca Y,particular(x) se alege sub forma unei dezvoltari în serie dupa functiile proprii Y(x), adica: si unde prin indicele k s-a identificat modul propriu de vibratie de ordinul k Tinând cont de faptul ca: ; dupa înlocuirea solutiei (7) în ecuatia de miscare (6) obtinem: Dezvoltând functiile trigonometrice si folosind definitiile anterioare ale marimilor n, ?, si ?, ultima ecuatie devine: Aceasta egalitate poate avea loc daca si numai daca valoarea coeficientului functiei cos(Ot) este nula, adica: ; Simplificând prin sin(Ot), se obtine: si, dupa înlocuirea definitiei functiei tg(?) si înmultirea cu (?A + ma): În vederea determinarii valorilor coeficientilor Ck vom înmulti ultima relatie cu Ym(x) si vom integra pe toata lungimea barei, L: si, folosind proprietatea de ortogonalitate a functiilor proprii: rezulta, în final: Pentru o forta perturbatoare constanta pe toata lungimea elementului combustibil, p0(x)=p0, coeficientii Ck devin: iar daca forta perturbatoare P0 este concentrata în sectiunea x = x*, valoarea corespunzatoare unitatii de lungime p0(x) se poate exprima cu ajutorul functiei simbolice Dirac, d(x): astfel încât coeficientii Ck devin: adica: Din formulele de mai sus rezulta ca marimile Ck sunt functii de pulsatia excitatoare O, valorile lor maxime fiind atinse pentru: Evident, în lipsa frecarilor interne si externe (µ = 0 si cv = 0) vom avea rezultatul imediat O = ?k Prin urmare, când pulsatia O a fortei perturbatoare este în aproprierea uneia din pulsatiile proprii ?k vom avea de-a face cu fenomenul de rezonanta, caracterizat prin cresterea semnificativa a amplitudinilor vibratiilor Fapt pentru care fiecarui coeficient Ck în parte i se asociaza factorul de amplificare: Forta excitatoare aleatoare Fara a intra în detalii privind calculele de analiza spectrala (metoda specifica utilizata în cazul semnalelor aleatoare), subliniem faptul ca daca prin F(O,t) se desemneaza forta excitatoare externa, atunci ecuatia ce descrie evolutia temporala a unui oscilator armonic cu amortizare vâscoasa (situatie valabila si cazului vibratiilor de încovoiere ale unei bare drepte si omogene) are forma generala: Aplicând acestei ecuatii transformata Fourier se obtine functia de raspuns în frecventa a sistemului analizat: unde, ca si în cazul analizat anterior, prin O se întelege o pulsatie caracteristica a fortei excitatoare, iar pin ?k pulsatia modului k de vibratie Daca sursa de excitatii aleatoare are o structura de zgomot alb (frecventele perturbatoare O se gasesc uniform distribuite în interiorul unei anumite benzi, limitata sa nu) si daca, în plus, S0 = constanta este densitatea spectrala a zgomotului alb, atunci bara (oscilatorul) va executa vibratii de încovoiere aleatoare tot de tip zgomot alb, densitatea spectrala a raspunsului sau fiind: Anulând derivata întâi a functiei S(O) în raport cu O avem: si rezulta ca functia S(O) prezinta maxime pentru pulsatiile excitatoare: Prin urmare, densitatea de putere a raspunsului barei va prezenta maxime (fenomen de rezonanta) pentru frecvente excitatoare O proportionale cu frecventele proprii ?k, valorile de extrem calculându-se cu formula: Alte elemente caracteristice raspunsului aleatoriu al barei, pe armonica k, vor fi: * raspunsul mediu (în timp) * raspunsul mediu patratic (în timp) * functia de autocorelatie unde, evident, prin Yk(x) s-a indicat functia fibrei deformate specifice armonicii de ordinul k VIBRATII INDUSE DE CURGEREA PARALELA Barele (tuburile, conductele etc ) aflate sub actiunea unui fluid curgând paralel cu acestea, prin exteriorul lor sau prin interiorul acestora, vor executa vibratii de încovoiere întretinute de forte putând avea caracter aleator (datorat turbulentei curgerii) sau pulsatoriu (datorat pulsatiilor de presiune generate de pompele de circulatie) În cazul vibratiilor induse de curgerea paralela a unui fluid se accepta existenta urmatoarelor forte ce actioneaza pe unitatea de lungime a barei studiate: * forta de antrenare transversala, dependenta de gradientul vitezei de curgere a fluidului; * forta de antrenare longitudinala, dependenta de modulul vitezei de curgere a fluidului; * fortele de presiune, normale pe suprafata laterala, rezultate din schimbarea fluxului momentan lateral al fluidului si considerate a fi independente de miscarea barei; * fortele de frecare externa, proportionale cu viteza de vibratie si opuse miscarii; * fortele de inertie, dependente de acceleratia vibratiilor si de masa elementului vibrator; * fortele elastice, dependente de rigiditatea barei si de derivata partiala de ordinul patru a deplasarii; * fortele de frecare interna, dependente de momentul de inertie si derivata partiala de ordinul doi a deplasarii; * fortele axiale de compresie sau întindere, dependente de caderea de presiune pe lungimea elementului studiat si modul de rezemare la capete Energia este creata de fortele de antrenare si disipata de celelalte sase descrise mai sus Prin urmare, primul pas în studierea raspunsului unei bare sub actiunea acestor forte consta în identificarea surselor de excitare În aceasta privinta se accepta ca mecanismele care participa cel mai frecvent si au un efect semnificativ în amplitudinea vibratiilor sunt cele care genereaza raspunsul aleator al sistemului elastic, ele putând fi împartite în patru categorii importante: * excitatia turbulenta; * excitatia fluido-elastica; * excitatia prin frecventa discreta; * excitatia datorata componentelor transversale ale curgerii paralele În functie de conditiile reale de functionare, sau de experimentare pe standuri de încercari, unele dintre mecanismele de excitare prezentate mai sus pot fi sursa dominanta de energie De cele mai multe ori însa este dificila identificarea excitatiei predominante datorita interactiunilor care pot avea loc între aceste mecanisme, descrise în cele ce urmeaza Excitatia turbulenta Fluctuatiile de presiune în curgerea turbulenta paralela pot cauza vibratii de încovoiere ale barelor si tuburilor Fluctuatiile apar în imediata vecinatate a peretilor rigizi sau sunt "transportate" spre acestia de catre perturbatiile curentului de fluid De obicei se accepta ca fluctuatiile de presiune care cauzeaza vibratii mecanice sunt generate: * în straturile turbulente de la capete, formându-se la suprafata rigida a barelor; * prin turbulenta rezultata din curgerea separata datorita eventualelor neuniformitati în canalul de Energiile asociate cu fluctuatiile de presiune turbulente sunt distribuite, în general într-o banda larga de frecventa si mai putin concentrate în jurul unei frecvente discrete Mecanismul excitatiilor turbulente consta în solicitari aleatoare dintre care bara selecteaza pentru raspunsul sau mai ales acea portiune de energie care este "închisa" în jurul unei frecvente proprii Vibratia nu trebuie sa cauzeze o crestere a fortelor turbulente efective, acestea fiind independente de miscarea vibratorie Excitatia fluido-elastica Vibratiile fluido-elastice sunt caracterizate printr-o viteza critica de curgere pâna la care amplitudinile de oscilatie sunt mici si dupa care acestea cresc foarte repede Sunt posibile doua mecanisme fluido-elastice: * un mecanism care depinde de pozitia sistemului vibrator (bara, tub etc ) si care este independent de timp; * un mecanism care depinde în special de timp În primul caz se considera ca exista în fiecare moment o situatie de curgere stationara Fortele care actioneaza aspra barei variaza atunci când pozitia sa relativa fata de fluid se schimba Datorita acestor schimbari, bara suporta forte mari de dragare, ce vor induce, astfel, vibratii Pentru mecanismul dependent de timp, fortele generate nu sunt numai functii de pozitia relativa a sistemului elastic fata de curgerea data, ci depind mai ales de existenta unei curgeri nestationare si de interactiunile temporare dintre bara vibratoare si câmpul de curgere Se presupune, de asemenea, ca deplasarea relativa a barei fata de fluid poate provoca si vibratii autoexcitate Excitatia prin frecvente discrete În general, fluctuatiile de presiune cauzate de separarea curgerilor prin valve, conducte si pompe sunt strâns legate de neuniformitatile în pasajul de curgere si sunt aleatoare, spectrul densitatii de putere al fluctuatiilor fiind distribuit într-o banda (mai larga sau mai îngusta de frecvente Totusi, exista si cazuri în care fluctuatiile de presiune pot fi periodice, cum ar fi, de exemplu, cele induse de functionarea pompelor de circulatie Daca frecventa unei componente a fluctuatiei este în imediata vecinatate a unei frecvente naturale a sistemului elastic, acesta va prezenta (executa) vibratii fortate de amplitudine mare (rezonanta) Nu toate frecventele discrete ale pulsatiilor aleatoare de presiune generate în sistemul bara /fluid vor excita vibratii de încovoiere De exemplu, unda de frecventa discreta si presiune aleatoare care se propaga în lungul unei bare poate cauza unde stationare axiale în sistem, dar nu poate cauza, în câmpul de curgere, gradienti transversali de presiune cantitativ apreciabili, astfel încât fortele transversale, actionând pe bara respectiva, tind sa se anuleze Aceste pulsatii vor trece în lungul barei ca o unda unidimensionala, viteza undei fiind egala cu viteza sunetului în fluid În fiecare moment, presiunea este egal distribuita în jurul suprafetei barei si nu contribuie semnificativ la vibratiile acesteia Excitatia componentei transversale de curgere Pentru barele aflate în curgere paralela si prezentând oscilatii de amplitudini relativ mici se accepta ca vibratiile sunt excitate în primul rând de componenta transversala a curgerii si de orice alta deviere de la curgerea stationara Aceste componente pot cauza vibratii ale barei, la fel ca si excitatiile de frecvente discrete, datorita raspândirii de vârtejuri de catre bara Prin urmare, cele doua mecanisme se pot suprapune, efectele lor însumându-se, fara a se putea specifica existenta predominanta a unuia sau altuia dintre ele Amortizarea vibratiilor de încovoiere Asa cum s-a aratat la începutul prezentului capitol, o parte din energia cedata de toate tipurile de mecanisme de excitare existente la un moment dat în curgerea paralela va fi disipata cu precadere prin amortizare Amortizarea totala a barei (tubului, conductei etc ) este egala cu suma dintre amortizarea hidrodinamica datorata fortelor de frecare externa, si amortizarea mecanica datorata fortelor de frecare interna Functie de situatiile reale experimentale sau de functionare, unul sau celalalt dintre aceste doua mecanisme poate fi predominant Amortizarea hidrodinamica Când o bara vibreaza într-un fluid, energia cinetica este comunicata fluidului si creste masa efectiva de oscilatie a barei Aceasta crestere, numita masa aditionala, este dependenta de geometria barei În cazul unui fascicul de bare, masa aditionala este suma maselor aditionale ale barelor componente Importanta practica a masei aditionale este aceea ca ea conduce la micsorarea valorilor frecventelor naturale Amortizarea hidrodinamica poate fi descompusa în: * amortizarea rezultata din vibratia barei în fluid stationar; * amortizarea rezultata din viteza de curgere axiala a fluidului de racire O tratare riguroasa a amortizarii hidrodinamice impune analiza interactiunilor dintre conditiile de frontiera, datorate constructiei, si cele de suprafata, cerute de curgerea axiala La viteze de curgere mici, primul tip de conditii este singura sursa de amortizare hidrodinamica, iar la viteze de curgere mari, conditiile de suprafata cerute de curgerea axiala devin cauze predominante ale amortizarii Evident, pentru viteze intermediare de curgere, ambele tipuri de conditii vor contribui la amortizare, interactiunea dintre ele fiind complexa Amortizarea mecanica În timpul vibratiilor, energia este disipata sub o forma sau alta si amplitudinea nu poate fi mentinuta constanta fara un aport continuu de energie O descriere a fortelor de amortizare mecanica asociate cu disiparea de energie este dificila, acestea putând fi functii de deplasare, viteza, solicitari sau alti factori În sfârsit, o alta forma de disipare mecanica a energiei miscarii vibratorii este cea cunoscuta sub numele de amortizare interna, caracteristica oricarui sistem elastic Raspunsul barei elastice la vibratii induse S-a vazut ca o bara introdusa într-un fluid aflat în curgere paralela poate vibra datorita unei varietati de cauze Raspunsul barei poate fi static sau dinamic Raspunsul static presupune o deformare statica, indusa de urmatoarele forte: * forte de frânare; * forte centrifuge; * forte date de presiunea statica a fluidului; * instabilitati statice Raspunsul dinamic presupune: * vibratii fortate; * vibratii autoexcitate si * rezonanta parametrica Vibratia fortata: la viteze de curgere subcritice, bara este supusa la excitatii de curgere Daca forta de excitatie este periodica, cum ar fi revarsarea vârtejurilor, respectiv, periodicitatea indusa de functionarea pompelor de circulatie, pot apare oscilatii cu amplitudini mari, strâns legate de fenomenul de rezonanta Daca excitatia este data de fluctuatii de presiune aleatoare, oscilatiile au, în general, amplitudini mici Vibratii autoexcitate: sunt vibratii cu amplitudine relativ mare, care au fost observate experimental, pentru viteze de curgere mari Rezonanta parametrica: curgerea pulsatorie si curgerea bifazica pot induce rezonante parametrice si combinatorii Vibratii de încovoiere în fluid stationar În cazul unei anumite structuri imersate într-un fluid, fluidul înconjurator trebuie sa fie deplasat pentru a usura miscarea Acest efect este luat în calcul prin intermediul conceptului de masa aditionala În mod ideal, masa aditionala poate fi determinata din ecuatia tridimensionala Navier-Stokes, dar un astfel de calcul este dificil chiar si în cazuri simple Pentru vibratii cu amplitudini mici, masa aditionala poate fi calculata cu ajutorul teoriei potentialului de curgere, teorie care are, totusi, unele limite Pentru un fluid vâscos, forta fluidului poate fi separata în doua componente: una în faza cu acceleratia (efectul de masa aditionala) si alta opusa miscarii (efectul de amortizare) În cele mai practice situatii, efectul vâscozitatii fluidului asupra masei aditionale este mic Însa vâscozitatea fluidului contribuie semnificativ la procesul de frânare, mai ales în cazurile în care amortizarea asociata cu vâscozitatea fluidului este foarte mare Atât marimea masei aditionale cât si amortizarea depind de amplitudinea vibratiilor Efectul compresibilitatii fluidului asupra miscarii structurii este similar cu cel al vâscozitatii deoarece, în cazul vibratiei unei bare singulare într-un volum infinit de fluid, masa aditionala depinde de lungimea de unda a perturbatiei ce se deplaseaza prin lichid La fel, amortizarea radiativa poate deveni semnificativa Cu alte cuvinte, compresibilitatea fluidului afecteaza modurile structurii datorita cuplajului fluid-bara elastica, fenomen asemanator cuplajului acustico-elastic Probleme mai complexe le ridica vibratiile cuplate ale elementelor multiplu structurale, cum este cazul unui fascicul de bare Datorita efectului de cuplaj al fluidului, miscarea unei bare într-un grup de mai multe bare va avea influente asupra celorlalte bare si fiecare bara va raspunde în cadrul grupului si mai putin individual Când un grup de cilindri vibreaza într-un fluid (vibratii cuplate), exista o relatie bine determinata pentru faza dintre cilindri Daca într-un grup de cilindri numai unul singur oscileaza si ceilalti ramân în repaus, miscarea primului se numeste necuplata Caracteristicile generale dinamice ale vibratiilor cuplate difera de cele pentru un singur cilindru Sa consideram, de exemplu, un grup de cilindri identici, cu o frecventa naturala în vid egala cu f0 si masa m Exista trei cazuri de vibratii care intereseaza: * Bara singulara în fluid infinit: * Vibratii necuplate pentru mai multe bare aflate în fluid: * Vibratii cuplate pentru mai multe bare aflate în fluid: unde: mf = masa de fluid deplasata; ma = masa aditionala; mp= valorile proprii ale matricei masei aditionale Vibratii induse de curgerea paralela Problema miscarii structurilor mecanice supuse unei curgeri externe sau/si continând un fluid în curgere este cunoscuta sub numele de vibratii induse de curgere În legatura cu orientarea curgerii fluidului fata de o anumita axa principala a structurii s-a convenit sa se clasifice problema vibratiilor induse în: * cazul curgerii paralele; * cazul curgerii transversale Fluidul aflat în curgere printr-un sistem este o sursa de energie care poate induce vibratii si instabilitati Fara curgere, sistemul se comporta ca o structura elastica în vid, dar frecventa este redusa datorita efectului masei aditionale În prezenta curgerii, sistemele conservative si cele neconservative raspund în mod diferit la cresterea vitezei de curgere În sistemele conservative frecventa creste cu descresterea vitezei de curgere Pentru o anumita valoare a vitezei de curgere, prima frecventa naturala tinde catre zero si sistemul îsi pierde stabilitatea, aparând fenomenul de flambaj Crescând, în continuare, viteza de curgere, alte frecvente naturale pot sa se anuleze Aceste viteze de curgere sunt asa-numitele viteze de curgere critice pentru modurile joase În sistemele neconservative frecventa este, în general, complexa La viteze de curgere mici sistemul executa oscilatii amortizate în toate modurile Alte efecte ale curgerii fluidului asupra evolutiei vibratiilor sunt: * reducerea frecventelor proprii; * contributii la amortizarea miscarii La anumite viteze de curgere mari, sistemul devine supra-amortizat Crescând, în continuare, viteza de curgere, sistemul poate prezenta instabilitati de tipul flambaj si/sau batai (zvâcniri) Un exemplu clasic de sistem neconservativ este cazul unui tub prin care curge fluid, tub care are capatul din amonte fix, iar cel din aval liber (grinda în consola) Daca si capatul din aval este fixat, sistemul poate fi considerat conservativ Pentru un fluid în curgere, ecuatiile de miscare contin termeni cu derivate mixte Din punct de vedere fizic, acestea sunt asociate cu fortele Coriolis Datorita prezentei lor sistemul este anizotropic, adica viteza de faza si viteza de grup sunt diferite în directia amonte fata de directia aval S-au corelat cu acceleratia Coriolis si alte fenomene interesante, cum ar fi * nu exista moduri normale în sisteme finite deoarece nu exista unde stationare; * pentru anumite valori ale frecventelor exista unde întârziate; * frecventa de taiere a unui mod de vibrare este coborâta de fluidul în miscare si nu corespunde unui numar de unda zero; * sistemul poate fi stabilizat sau destabilizat de catre fortele Coriolis Si efectul presiunii statice a fluidului poate fi semnificativ, prin aparitia unei deformatii statice În bare va exista o forta laterala (transversala) a carei actiune este echivalenta unei compresii axiale Astfel, crescând presiunea statica, se tinde sa se reduca frecventa de oscilatie si, prin presurizare, se poate induce flambajul Într-un sistem constituit din mai multe structuri, cum ar fi, de exemplu, un fascicul de elemente combustibile, apar interactiuni între acestea Sunt importante doua efecte: * datorita cresterii masei aditionale, frecventele si vitezele critice de curgere sunt mai mici decât cele ale unei singure bare supuse unei curgeri paralele; * frecventele si vitezele critice de curgere pentru sistemul cuplat sunt mai mici decât cele pentru un tub cilindric prin care trece un fluid (tubul de presiune) Vibratii subcritice si instabilitati Pentru întelegerea vibratiilor subcritice s-au sugerat trei mecanisme: * vibratii autoexcitate; * vibratii fortate; * rezonanta parametrica Cercetarile efectuate au condus însa la concluzia ca modelul vibratiilor fortate este cel mai potrivit În general, este acceptat faptul ca o structura supusa unei curgeri externe suporta zgomotul de curgere, care devine primul mecanism (predominant) al excitatiei fortate Excitatia parametrica apare numai într-o curgere stationara monofazica, iar vibratia autoexcitata poate apare si într-o curgere bifazica, dar efectele lor sunt mult mai reduse în cazul vibratiilor subcritice Fenomene foarte interesante, asociate cu sistemele neconservative, au condus la o serie de cercetari ale sistemului neconservativ tub-fluid Un subiect controversat asociat cu un astfel de sistem este cel cunoscut ca efectul de destabilizare al fortelor dependente de viteza Viteza de destabilizare obtinuta din ecuatia de miscare în absenta fortelor dependente de viteza este mai mica decât viteza de destabilizare obtinuta prin includerea în calcule si a fortelor dependente de viteza Este un paradox, consecinta a aplicarii necritice a metodei micilor vibratii De aceea s-a introdus conceptul de cuazi-stabilitate În cazurile când toate frecventele sunt reale, se arata ca exista o parte considerabila a regiunii cuazi-stabile, îndepartata de regiunea de instabilitate, în care actioneaza fortele dependente de viteza Doua proprietati interesante în sistemele neconservative structura-fluid sunt: * efectul de stabilizare al conditiilor de rezemare; * efectul de batai Alaturi de flambaj si batai mai exista si alte forme de instabilitate în curgerea paralela Astfel, în curgerile nestationare si curgerile bifazice, excitatiile parametrice, datorate schimbarilor periodice ale vitezei de curgere si ale masei, pot induce rezonante parametrice si rezonante combinate Acest lucru s-a observat în curgerea pulsata, pentru sisteme continue sau discrete, care sunt stabile pâna la o anumita viteza de curgere si care îsi pot pierde stabilitatea prin rezonante parametrice si combinate Mai mult, daca viteza de curgere creste, surplusul de energie cinetica cauzeaza si el o stare de instabilitate Am vazut, deci, ca vibratiile barelor se pot exprima ca moduri normale de încovoiere, ale caror surse (în general aleatoare, dar nu numai) sunt, în special, mecanismele de inducere prin curgerea paralela stabilizata a fluidului de racire O analiza riguroasa a raspunsului unei bare la astfel de excitatii impune utilizarea analizei Fourier si a conceptelor statistice Bara vibreaza având un numar infinit de frecvente naturale si moduri asociate de vibratie, fapt pentru care, în analizele statistice, se cere folosirea a cât mai multor moduri Totusi, deoarece modul fundamental este predominant, se pot accepta simplificarile date de folosirea în calcule numai a acestui mod În acelasi timp, amplitudinea de vibratie a barei studiate poate fi descompusa în doua componente pe doua directii perpendiculare Astfel, pentru surse de excitare uniform distribuite pe directie circumferentiala, o bara cu conditii de rezemare identice la ambele capete va vibra într-o maniera statistica similara în orice plan De aceea este suficient sa se considere / analizeze doar componenta pe o singura directie Altfel spus, componenta fortei aleatoare pe o directie va fi folosita la calcularea amplitudinii vibratiei numai în directia respectiva Acest lucru permite identificarea barei cu un oscilator cu un singur grad de libertate Rezulta, prin urmare, necesitatea caracterizarii cât mai exacte a câmpului de presiune din jurul barei, ca si identificarea, pe cât posibil, a mecanismelor de inducere a vibratiilor de încovoiere pentru diferite cazuri concrete de functionare (sau de experimentare) Concordanta dintre datele experimentale si cele obtinute dintr-un eventual model teoretic al situatiei reale studiate va depinde si de îndeplinirea acestor ultime doua cerinte O generalizare a ecuatiei de miscare Pentru exemplificarea celor prezentate mai sus, dam în continuare ecuatia de miscare pentru o bara supusa excitatiilor induse de curgerea paralela a unui fluid oarecare Având viteza medie V, curgerea paralela a fluidului actioneaza asupra unui element de lungime ?x cu urmatoarele forte suplimentare: * forta laterala datorata schimbarii momentului lateral al curgerii: * forta de antrenare transversala: * forta de antrenare longitudinala: unde: ?f = densitatea fluidului; D = diametrul exterior al elementului combustibil; CM = coeficient de masa aditionala; CD = coeficient de antrenare laterala; Cf = coeficient de antrenare longitudinala În ce priveste forta de tensionare N, aceasta este suma dintre tensiunea de montaj, N0, si o forta de tensiune datorata frecarii longitudinale cu fluidul aflat în cuegere paralela: unde prin C'T s-a notat coeficientul de antrenare din partea fluidului În aceste conditii vom avea: cu ? = 1 daca capatul din aval este simplu sprijinit si ? = 0 daca în capatul din aval avem reazem elastic sau capat liber Cu aceste noi definitii, ecuatia vibratiilor de încovoiere ale barei aflate în curgere paralela devine: unde q(x,t) reprezinta forta exercitata asupra elementului combustibil din partea turbulentei curgerii Fara a insista asupra modului de rezolvare a ecuatiei de mai sus, vom prezenta doar cele mai importante concluzii calitative ce se pot desprinde dintr-o asemenea analiza: * functie de conditiile de rezemare la capete, frecventa fundamentala poate sa creasca sau sa descreasca cu viteza de curgere; pentru o bara cu capetele fixate, frecventa fundamentala descreste cu cresterea vitezei de curgere a fluidului, iar pentru o bara cu capetele libere, frecventa fundamentala creste odata cu cresterea vitezei fluidului; * masa aditionala creste considerabil daca bara are ambele capete fixate; mai mult, în cazul unui fascicul de bare, trebuie luat în considerare si efectul barelor adiacente; * amortizarea vibratiilor este amplificata de cresterea vitezei de curgere a fluidului, datorita fortelor normale de antrenare si fortelor Coriolis; * amplitudinea vibratiilor creste cu viteza de curgere si poate fi exprimata printr-o functie exponentiala, exponentul depinzând de valoarea fortelor excitatoare, de amortizare, de diametrul hidraulic si de conditiile de rezemare la capete; * crescând rigiditatea reazemelor de la capetele barei se poate reduce valoarea elongatiilor vibratiilor de încovoiere CONCLUZII Pastrând simetria dimensiunilor primului si ultimului capitol, subliniem doar faptul ca, în general, studiile teoretice de Rezistenta materialelor se bazeaza pe modelul idealizat de bara Bernoulli-Euler (lungimea " aria sectiunii transversale, oscilatii mici, etc ) Luând în calcul si distributia parabolica a fortelor de forfecare în sectiunea transversala, pentru vibratiile libere ale asa-numitei bare Timoshenko cu rigiditate EI constanta se deduce ecuatia de miscare: unde R = I / A si K este un coeficient de forfecare prin care se introduce distributia parabolica De remarcat prezenta derivatei partiale de ordinul 4 în raport cu timpul BIBLIOGRAFIE 1 ANDO Y Three-dimensional theory on the natural vibration of rods ; Nuclear Engineering and Design Vol 15 (1971) 2 BELYTSCHKO T Methods and programs for analysis of fluid-structure systems ; Nuclear Engineering and Design, Vol 42 (1977) 3 BELYTSCHKO T Fluid-structure interactions in light water reactor systems ; Nuclear Engineering and Design, Vol 60 (1980) 4 BUZDUGAN GH Rezistenta materialelor ; Editura Tehnica, Bucuresti, 1980 5 BUZDUGAN GH , FETCU L , RADES M Vibratii mecanice ; Editura Didactica si Pedagogica, Bucuresti, 1979 6 CAIUS I Mecanica teoretica ; Editura Didactica si Pedagogica , Bucuresti, 1980 7 CARLUCCI L N Damping and hydrodynamic mass of a cylinder in two-phase flow ; Atomic Energy of Canada Limited-7315, 1980 8 CHEN S S Parallel-flow-induced vibration of fuel rods ; Nuclear Engineering and Design, Vol 18 (1972) 9 DEUTSCH I Rezistenta Materialelor ; Editura Didactica si Pedagogica, Bucuresti, 1979 10 DOCA C On the Vibration / Wear Correlations ; Journal of Materials Science and Technology, Vol 8-No 1/2000, pag 109-114, Coral Press Scientific Publishing, Sofia 11 DOCA C Two Theoretical Models for "Flow Assisted Fueling" Changing Sequence in a CANDU Reactor ; Journal of Materials Science and Technology, Vol 9-No 1/2001, pag 3-11, Coral Press Scientific Publishing, Sofia 12 DOCA C , DOCA L Extrapolation of Monotonous Functions using the Linear Correlation Coefficient ; - Scientific Bulletin of "Politehnica", University of Bucharest, Volume 61, 1999, Number 1-2, pag 177-184 13 DOCA C , DOCA L A Model of the Laminar Flow with Phase Exchange for Real Liquids Quantitative Measurement of the H2O Losses through the Imperfect Tightness in a High Pressure and Temperature Installation ; Scientific Bulletin of "Politehnica" University of Bucharest, Volume 61, 1999, Number 1-2, pag 185-194 14 DOCA C, DOCA L , BÎRNEANU I On the Fretting Wear in CANDU Pressure Tube ; Scientific Bulletin of "Politehnica" University of Bucharest, Serie A, Vol 62-No 2/2000, pag 53-61 15 DOCA C , PANESCU M On the Thermal Conductivity of Sodium Vapour ; Studia Univ Babes - Bolyai Cluj- Napoca, Physica, XXXIV, 1, 1990, pag 53-61 16 ELVIS P D Flow Induced Vibrations ; Von Fostrand Reinhals Co 1977 17 FERTIS D G Dynamics and Vibrations of Structures ; Jhon Wiley & Soons, New York, 1973 18 GORDIS K Outline of dynamic analysis for piping systems ; Nuclear Engineering and Design, Vol 52 (1979) 19 KERLIN T W Frequency response testing in nuclear reactor ; Academic Press, New York, 1974 20 KUZELKA V Vibrational response of reactor components to random forces ; Nuclear Engineering and Design, Vol 72 (1982) 21 MANGERON D , IRIMICIUC N Mecanica rigidelor cu aplicatii în inginerie Mecanica vibratiilor sistemelor de rigide ; Editura Tehnuica, Bucuresti, 1981 22 MASON V A Dynamics of nuclear fuel assemblies in vertical flow channels ; Atomic Energy of Canada Limited-5976, 1978 23 PAIDOUSSIS M P Dynamics of flexible cylinders in axial flow ; Atomic Energy of Canada Limited-6387, 1979 24 PETTGREW M J Dynamics and stability of flexible cylinders subjected to axial flow; Atomic Energy of Canada Limited-5502, 1976 25 PETTIGREW M J Flow induced vibration of nuclear power station components ; Atomic Energy of Canada Limited-5852, 1977 26 PETTGREW M J Flow- induced vibration in nuclear power station components ; Atomic Energy of Canada Limited-7291, 1981 27 POSEA N Rezistenta materialelor ; Editura Didactica si Pedagigica, Bucuresti, 1979 28 RADES M Metode dinamice pentru identificarea sistemelor mecanice ; Editura Acaddemiei RSR, 1979 29 ROSINGER H E , RITCHIE I G On Timoshenko'correction for shear in vibrating isotropic beams, Atomic Energy of Canada Limited- 57644, 1977 30 ROSCULET M Analiza matematica ; Editura Didactica si Pedagogica Bucuresti, 31 TAKEUCHI K Evaluation of hydrodynamic mass methods; Nuclear Engineering and Design, Vol 70 (1982) ANEXE ANEXA A - Functii hiperbolice Definitii si identitati: ; ; Formulele lui Euler: Relatii dintre functiile hiperbolice si cele trigonometrice: ANEXA B - Notatii si unitati de masura Semnfificatie Simbol Unitate de masura Axe de coordonate, cu origine în centrul geometric al sectiunii transversale de arie A x, y, z m Moment de încovoiere M Nm Forta normala N N Forta de forfecare pe directia i (= y, z) Ti N Sarcina pe unitatea de lungime q(x) N/m Forta externa concentrata F N Aria sectiunii transversale A m2 Lungime L m Deplasarea / deformarea pe directia x u(x) m Deplasarea / deformarea pe directia y v(x) m Sageata de încovoiere pe directia z w(x) m Moment de inertie (moment de ordinul II al sectiunii transversale) Iy, Iz m4 Modul de elasticitate (modulul lui Young) E N/m2 Coeficientul lui Poisson µ - Modul de forfecare G N/m2 ANEXA C - Momente de inertie Daca x, y si z sunt coordonate carteziene, atunci momentele de inertie (momentele de ordinul doi) ale unei sectiuni de arie A, perpendiculare pe axa x se definesc: ; ; Tabel 1 - Exemple de formule de calcul al momentului de inertie Sectiune transversala Momente de inertie ; ; ; ; ; ANEXA D - Rezemarea la capete Plecând de la definitiile reazemelor de capat "consacrate" (fara forte concentrate si momente încovoietoare în abscisele x0 = 0 si xL = L), fizic sunt posibil urmatoarele 37 de combinatii pentru o bara orizontala: Tabel 2 - Exemple de rezemarea la capete Tip de rezemare / Simbol Ecuatii de caracterizare Capat liber - încastrare rigida Capat liber - încastrare elastica Capat liber - rezemare generala Reazem alunecator - încastrare rigida Reazem alunecator - încastrare elastica Reazem alunecator - rezemare generala Reazem simplu rigid - reazem simplu rigid Reazem simplu rigid - încastrare rigida Reazem simplu rigid - reazem simplu elastic Reazem simplu rigid - încastrare elastica Reazem simplu rigid - rezemare generala Încastrare rigida - capat liber Încastrare rigida - reazem alunecator Încastrare rigida - reazem simplu rigid Încastrare rigida - încastrare rigida Încastrare rigida - reazem simplu elastic Încastrare rigida - încastrare elastica Încastrare rigida - rezemare generala Reazem simplu elastic - reazem simplu rigid Reazem simplu elastic - încastrare rigida Reazem simplu elastic - reazem simplu elastic Reazem simplu elastic - încastrare elastica Reazem simplu elastic - rezemare generala Încastrare elastica - capat liber Încastrare elastica - reazem alunecator Încastrare elastica - reazem simplu rigid Încastrare elastica - încastrare rigida Încastrare elastica - reazem simplu elastic Încastrare elastica - încastrare elastica Încastrare elastica - rezemare generala Rezemare generala - capat liber Rezemare generala - reazem alunecator Rezemare generala - reazem simplu rigid Rezemare generala - încastrare rigida Rezemare generala - reazem simplu elastic Rezemare generala - încastrare elastica Rezemare generala - rezemare generala ANEXA E - Calculul constantelor de integrare pentru o bara fara reazeme intermediare, în cel mai simplu caz de solicitare: fara pat elastic, fara forta axiala, moment de inertie constant, sarcina distribuita constanta Asa cum s-a putut vedea în capitolul ECUATII, ALGORITMI SI SOLUTII, prin rezolvarea oricareia dintre cele 36 de ecuatii de încovoiere analizate se obtine câte un set corespunzator de patru constante de integrare: C0, C1, C2 si C3, necunoscute Acelasi set de patru constante se obtine si în cazul vibratiilor de încovoiere Aceste necunoscute se evalueaza plecând de la ecuatiile de caracterizare ale reazemelor de capat si, daca ne restrângem doar la cele 37 de situatii din Tabelul 2, ar rezulta cel putin 36 × 37 = 1332 de descrieri complete privind încovoiere a unei bare studiate S-a subliniat sintagma cel putin deoarecea acestor tipuri de probleme li se mai poate adauga înca o multime de alte combinatii înglobând, pe de o parte reazeme "consacrate" si, pe de alta parte, actiunea externa în capatele barei fie a unei forte concentrate, fie a unui moment încovoietor, fie a amândurora Din acest noian de posibilitati, vom prezenta în continuare, ca exemplu de frumusete (si complexitate) a calcului analitic, doar cele mai simple situatii ce se încadreaza cazului de solicitare identificat în titlul prexentei ANEXE Astfel, constantele de integrare din solutia: respectiv din derivatele: se calculeaza impunând acestor functii sa satisfaca conditiilor corespunzatoare naturii reazemelor de capat din abscisele x0 = 0 si xL = L, dupa cum urmeaza: * Capat liber - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Capat liber - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Capat liber - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem alunecator - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem alunecator - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem alunecator - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu rigid - reazem simplu rigid Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu rigid - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu rigid - reazem simplu elastic Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu rigid - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu rigid - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - capat liber Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - reazem alunecator Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - reazem simplu rigid Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - reazem simplu elastic Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare rigida - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu elastic - reazem simplu rigid Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu elastic - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu elastic - reazem simplu elastic Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu elastic - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Reazem simplu elastic - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - capat liber Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - reazem alunecator Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - reazem simplu rigid Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - reazem simplu elastic Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Încastrare elastica - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - capat liber Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - reazem alunecator Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - reazem simplu rigid Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si: * Rezemare generala - încastrare rigida Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - reazem simplu elastic Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - încastrare elastica Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si * Rezemare generala - rezemare generala Rezolvând în necunoscutele C0, C1, C2 si C3 sistemul format din cele patru ecuatii de caracterizare a reazemelor de capat: se obtin solutiile: si ANEXA F - LICENTIEREA GNU FDL Aceasta este o traducere neoficiala în limba româna a Licentei GNU pentru Documentatie Libera (GNU Free Documentation Licence) Textul original al licentei, în limba engleza, se gaseste si poate fi descarcat de la: http://www gnu org/copyleft/fdl html Licenta GNU pentru Documentatie Libera (GNU Free Documentation Licence) Versiunea 1 2 din noiembrie 2002 Copyright (C) 2000, 2001, 2002 Free Software Foundation, Inc 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA Oricine poate copia si distribui copii identice ale acestui document, dar modificarea lui nu este permisa 0 PREAMBUL Scopul acestei Licente este de a conferi unui set de instructiuni, manual scolar sau altui document folositor "libertate", înteleasa în sensul urmator: asigura tuturor libertatea de a copia si redistribui textul, cu sau fara modificari, în scopuri comerciale si necomerciale Ca scop secundar, aceasta Licenta rezerva pentru autor si editor dreptul de a fi creditati pentru munca lor si de a nu fi responsabili pentru modificarile efectuate de altii Aceasta Licenta confera un fel de "copyleft", ceea ce înseamna ca lucrarile derivate trebuie sa fie si ele libere în sensul de mai sus Aceasta Licenta este inspirata de Licenta Publica Generala GNU (GNU General Public Licence, GNU GPL), care este o licenta similara conceputa pentru a acoperi softul liber Aceasta Licenta a fost scrisa pentru a acoperi manuale pentru soft liber, pentru ca softul liber necesita documentatie libera: un program trebuie însotit de manuale care ofera aceeasi libertate în folosire ca si softul Aceasta Licenta nu este limitata la manuale pentru soft si poate fi folosita pentru a acoperi orice lucrare, indiferent de subiect sau de modul de publicare Aceasta Licenta este recomandata în principal pentru lucrari care servesc drept referinta sau au fost scrise în scop de instruire 1 APLICABILITATE SI DEFINITII Aceasta Licenta se aplica oricarui manual sau lucrari, în orice mediu, care contine o nota inclusa de catre detinatorul dreptului de autor ce permite distributia sub acoperirea acestei Licente Aceasta nota confera dreptul universal (world-wide), fara indemnizatie si nelimitat ca durata de a folosi aceasta lucrare în conditiile descrise de aceasta Licenta Termenul "Documentul" folosit mai jos se refera la manualul sau lucrarea acoperita de Licenta Orice membru al publicului este un beneficiar al acestei Licente si va fi desemnat prin termenul "Dvs " sau prin folosirea persoanei a doua Se considera în mod automat ca ati acceptat termenii acestei Licente daca copiati, modificati sau distribuiti Documentul într-un mod ce necesita permisiunea autorului în conformitate cu legea drepturilor de autor O "Versiune Modificata" a Documentului este orice lucrare continând Documentul sau o portiune din Document, copiata identic sau cu modificari si / sau tradusa într-o alta limba O "Sectiune Secundara" este o anexa cu titlu, sau o sectiune mentionata în cuprins care are ca scop exclusiv descrierea relatiei editorilor sau a autorilor Documentului cu subiectul Documentului (sau cu subiecte legate de acesta) si care nu contine subiecte incluse în mod direct în subiectul Documentului (Asadar, daca Documentul este în parte manual de matematica, o Sectiune Secundara nu poate contine explicatii matematice ) Relatia poate fi o conexiune istorica cu subiectul sau cu problemele înrudite cu subiectul, sau puncte de vedere legale, comerciale, filozofice, etice sau politice legate de acesta "Sectiunile Invariante" sunt anumite Sectiuni Secundare ale caror titluri sunt specificate ca fiind titluri de Sectiuni Invariante din Document în nota ce permite distributia Documentului sub acoperirea acestei Licente Daca o sectiune nu este conforma cu definitia de mai sus a unei Sectiuni Secundare ea nu poate fi desemnata drept Sectiune Invarianta Documentul poate sa nu contina nici o Sectiune Invarianta Daca Documentul nu specifica nici o Sectiune Invarianta se considera ca nu exista nici una "Textele De Coperta" sunt pasaje scurte de text care sunt listate ca Texte Pentru Coperta I (coperta din fata) si Texte Pentru Coperta IV (coperta din spate) în nota ce permite distributia Documentului sub acoperirea acestei Licente Un Text Pentru Coperta I poate avea cel mult 5 cuvinte, iar un Text Pentru Coperta IV poate avea cel mult 25 de cuvinte O copie "Transparenta" a Documentului este o copie în format electronic, reprezentata într-un format a carui specificatie este disponibila publicului, care este usor de modificat folosind un editor de text generic sau (pentru imagini compuse din pixeli) un editor grafic generic sau (pentru desene) un editor larg raspândit de grafica vectoriala, si care poate fi folosit de catre programe de formatare de text sau de catre programe de conversie în alte formate care pot fi folosite ca intrare de catre programe de formatare a textului O copie facuta într-un format de fisier Transparent dar care prin prezenta sau absenta anumitor elemente specifice formatului descurajeaza sau împiedica modificarile ulterioare nu este o copie Transparenta Un format grafic - o imagine - nu este un format Transparent daca este folosit pentru a reprezenta o cantitate substantiala de text O copie care nu este "Transparenta" este "Opaca" Exemple de formate compatibile cu copiile Transparente includ: text ASCII fara marcare, format de intrare Texinfo, format de intrare LaTeX, SGML si XML folosind un DTD public, HTML simplu si standard, fisiere PostScript si PDF modificabile Exemple de formate Transparente pentru imagine includ PNG, XCF si JPG Formatele Opace includ formate de text ce pot fi citite si editate doar de procesoare de text particulare (proprietary), SGML si XML pentru care DTD-ul si / sau uneltele de procesare nu sunt disponibile, HTML generat automat, documente PostScript si PDF produse de diverse procesoare de text exclusiv în scopul printarii / afisarii "Pagina de Titlu" înseamna, pentru o carte tiparita, pagina cu titlul si paginile urmatoare necesare pentru a tipari lizibil materialul care trebuie tiparit conform acestei Licente pe Pagina de Titlu Pentru lucrari care nu au o pagina cu titlu propriu-zisa "Pagina de Titlu" este textul aflat lânga principala aparitie a titlului lucrarii, precedând începutul corpului Documentului O sectiune "Numita XYZ" este o sectiune din Document al carei titlu este fie XYZ sau contine XYZ în paranteze dupa textul care traduce XYZ în alta limba (Aici XYZ înlocuieste nume specifice ce vor fi mentionate mai jos, ca de exemplu "Multumiri", "Dedicatii", "Giruri" (Endorsement) si "Istorie" ) A "Pastra Titlul" unei astfel de sectiuni atunci când modificati Documentul înseamna ca aceasta ramâne "Numita XYZ" conform acestei definitii Documentul poate include Limitari de Responsabilitate (Warranty Disclaimers) atasate notificarii care afirma ca aceasta Licenta se aplica Documentului Aceste Limitari de Responsabilitate se considera a fi incluse pentru referinta în aceasta Licenta: orice alte implicatii pe care aceste Limitari de Responsabilitate le-ar putea avea sunt nule si nu au nici un efect asupra întelesului acestei Licente 2 COPII IDENTICE Puteti copia si distribui Documentul pe orice mediu, comercial sau necomercial, atâta timp cât aceasta Licenta, notificarile de drepturi de autor si notificarea de licenta care spune ca aceasta Licenta se aplica acestui Document sunt reproduse în toate copiile, si atâta timp cât nu adaugati nici un fel de alta conditie în afara de cele prezente în aceasta Licenta Nu aveti dreptul sa luati masuri tehnice de a obstructiona sau controla citirea sau recopierea copiilor pe care le faceti sau le distribuiti Aveti totusi dreptul sa acceptati compensatii în schimbul copiilor Daca distribuiti un numar suficient de mare de copii trebuie sa respectati si conditiile din sectiunea 3 Aveti de asemenea dreptul sa împrumutati copii în aceleasi conditii ca cele de mai sus, si aveti dreptul sa afisati copii 3 COPIEREA ÎN CANTITATI MARI Daca publicati copii tiparite (sau copii în medii care folosesc de obicei coperte tiparite) ale Documentului, în numar mai mare de 100 si daca notificarea de licenta a Documentului cere Texte de Coperta, trebuie sa includeti copiile pe coperte care sa contina, clar si lizibil, toate aceste Texte de Coperta: Textele Pentru Coperta I pe coperta I si Texte Pentru Coperta IV pe coperta IV Ambele coperte trebuie de asemenea sa va identifice în mod clar si lizibil ca editor al respectivelor copii Coperta I trebuie sa prezinte titlul în întregime, cu toate cuvintele din titlu la fel de vizibile si proeminente Puteti adauga alte materiale pe coperta în plus Copierea cu modificarile limitate la coperte, atâta timp cât satisfac aceste conditii, pot fi tratate în toate celelalte aspecte ca si copii identice Daca textele necesare pentru oricare dintre coperte sunt prea voluminoase pentru a încapea în mod lizibil, trebuie sa le includeti pe primele în ordinea originala (atâtea câte încap în mod rezonabil) pe coperta efectiva si sa continuati cu restul pe pagini adiacente Daca publicati sau distribuiti copii Opace ale documentului în numar mai mare de 100, trebuie ori sa includeti câte o copie Transparenta în format electronic împreuna cu fiecare copie Opaca, sau sa specificati în sau împreuna cu fiecare copie Opaca o locatie de retea electronica la care publicul general care foloseste reteaua sa aiba acces pentru a descarca, folosind un protocol standard public, copii complete Transparente ale documentului, fara adaugarea oricarui material aditional Daca folositi a doua optiune trebuie sa faceti demersuri rezonabil de prudente ca atunci când începeti distribuirea copiilor Opace sa va asigurati ca aceasta copie Transparenta va ramâne accesibila în acest fel la locatia respectiva timp de cel putin un an dupa distributia ultimei copii Opace (în mod direct sau prin agenti sau distribuitori) a acelei editii pentru public Se cere, dar nu în mod necesar, sa contactati autorii Documentului cu o perioada buna înainte de a distribui orice cantitate mare de copii, pentru a le da ocazia sa va puna la dispozitie o versiune actualizata a Documentului 4 MODIFICARI Puteti copia si distribui o Versiune Modificata a Documentului în conditiile sectiunilor 2 si 3 de mai sus, cu conditia de a acoperi Versiunea Modificata sub exact aceasta Licenta, cu Versiunea Modificata tinând locul Documentului, astfel licentiind distribuirea si modificarile Versiunii Modificate oricui intra în posesia unei copii ale acesteia În plus, trebuie sa faceti urmatoarele lucruri în Versiunea Modificata: A) Folositi în Pagina de Titlu (si pe coperte, daca exista) un titlu diferit de cel al Documentului, si de versiunile sale anterioare (care trebuie, daca exista, sa fie listate în sectiunea de Istorie a Documentului) Puteti folosi acelasi titlu ca o versiune anterioara daca editorul original al acelei copii va da permisiunea B) Listati pe Pagina de Titlu, ca autori, una sau mai multe dintre persoanele sau entitatile responsabile în calitate de autori pentru modificarile Versiunii Modificate, împreuna cu cel putin cinci dintre autorii principali ai Documentului (toti autorii principali, daca are mai putin de cinci), în afara de cazul ca acestia va elibereaza de aceasta obligatie C) Includeti pe Pagina de Titlu numele editorului Versiunii Modificate în calitate de editor D) Pastrati toate notificarile de drepturi de autor ale Documentului E) Adaugati o notificare de drepturi de autori relevanta pentru modificarile Dvs adiacent celorlalte notificari de drepturi de autor F) Includeti, imediat dupa notificarile de drepturi de autor, o notificare de licenta dând permisiune publica de a folosi Versiunea Modificata în conditiile acestei Licente, sub forma prezentata în Apendicele de mai jos G) Pastrati în acea notificare de licenta lista integrala a Sectiunilor Invariante si Textele de Coperta necesare date în notificarea de licenta a Documentului H) Includeti o copie nealterata a acestei Licente I) Pastrati sectiunea Numita "Istorie", Pastrati-i Titlul si adaugati-i un element care sa indice macar titlul, anul, noii autori si editorul Versiunii Modificate asa cum este dat pe Pagina de Titlu Daca nu exista o sectiune Numita "Istorie" în Document, creati una în care indicati titlul, anul, autorii si editorul Documentului asa cum este dat pe Pagina de Titlu al acestuia si apoi adaugati un element care sa descrie Versiunea Modificata asa cum a fost cerut în fraza precedenta J) Pastrati locatia de retea, daca exista, data în Document pentru acces public la o copie Transparenta a Documentului, cât si locatiile de retea date în Document pentru versiunile mai vechi pe care s-a bazat acesta Acestea pot fi incluse în sectiunea Numita "Istorie" Puteti omite locatia de retea a unei lucrari care a fost publicata cu cel putin patru ani înainte de Documentul în sine, sau daca editorul original al versiunii la care se refera va da permisiunea K) Pentru orice sectiune Numita "Multumiri" sau "Dedicatii" Pastrati Titlul sectiunii si pastrati în sectiunile respective toata substanta si tonul multumirilor si dedicatiilor fiecarui contribuitor L) Pastrati toate Sectiunile Invariante ale Documentului, nealterate ca text si ca titluri Numerotarea sectiunilor sau echivalentul numerotarii nu sunt considerate ca facând parte din titlurile sectiunilor M) Stergeti orice sectiune Numita "Giruri" O astfel de sectiune nu poate fi inclusa în Versiunea Modificata N) Nu modificati titlul nici unei sectiuni existente pentru a fi Numita "Giruri" sau pentru a intra în conflict cu vreo Sectiune Invarianta O) Pastrati toate Limitarile de Responsabilitate Daca Versiunea Modificata include sectiuni noi incluse în titlu sau anexe care se califica drept Sectiuni Secundare si nu contin material copiat din Document, aveti dreptul la alegerea Dvs sa numiti unele sau toate acestea ca fiind sectiuni invariante Pentru a face aceasta, adaugati-le titlurile la lista de Sectiuni Invariante în notificarea de licenta a Versiunii Modificate Aceste titluri trebuie sa fie distincte fata de toate celelalte titlurile de sectiune Puteti adauga o sectiune Numita "Giruri" doar daca aceasta contine numai girurile a diverse entitati asupra Versiunii Modificate - de exemplu recenzii sau faptul ca textul a fost aprobat de o organizatie ca fiind o definitie autoritara a unui standard Puteti adauga un pasaj de cel mult cinci cuvinte ca Text Pentru Coperta I si un pasaj de cel mult 25 de cuvinte ca Text Pentru Coperta IV la sfârsitul Textelor De Coperta în Versiunea Modificata Numai un singur pasaj poate fi adaugat la Textul Pentru Coperta I si unul la Textul Pentru Coperta IV de catre (sau prin aranjament cu) orice entitate Daca Documentul contine deja texte de coperta pentru coperta respectiva, adaugat în prealabil de Dvs sau prin aranjament cu aceeasi entitate în numele careia actionati, atunci nu puteti adauga un altul, însa puteti sa-l înlocuiti pe cel vechi numai cu permisiunea explicita a editorului anterior care l-a adaugat pe cel vechi Autorul (autorii) si editorul (editorii) Documentului nu va dau prin aceasta Licenta permisiunea de a le folosi numele pentru publicitate sau pentru a pretinde sau implica vreo girare a oricarei Versiuni Modificate 5 COMBINAREA DOCUMENTELOR Puteti combina Documentul cu alte documente acoperite de aceasta Licenta sub termenii definiti în sectiunea 4 de mai sus pentru versiuni modificate, cu conditia sa includeti în versiunea combinata toate Sectiunile Invariante ale tuturor documentelor originale, nemodificate, si sa le listati pe toate ca Sectiuni Invariante ale versiunii combinate în notificarea de licenta, cât si sa pastrati toate Limitarile de Responsabilitate Versiunea modificata nu trebuie sa contina decât o singura copie a acestei Licente, iar duplicatele identice ale Sectiunilor Invariante pot fi înlocuite cu o singura copie Daca exista Sectiuni Invariante cu nume identice si continut diferit, schimbati-le numele adaugând la sfârsitul titlului, în paranteza, ori numele autorului sau al editorului original al acelei sectiuni daca acesta este cunoscut, ori un numar unic Faceti aceleasi modificari respective titlurilor sectiunilor în lista de Sectiuni Invariante din notificarea de licenta a versiunii combinate În versiunea combinata trebuie sa combinati si toate sectiunile Numite "Istorie" din diversele documente originale, creând o sectiune unica Numita "Istorie"; la fel trebuie sa combinati si toate sectiunile Numite "Multumiri" cât si cele Numite "Dedicatii" Trebuie sa stergeti toate sectiunile Numite "Giruri" 6 COLECTII DE DOCUMENTE Puteti crea o colectie formata din Document si alte documente acoperite de aceasta Licenta si sa înlocuiti copiile individuale ale acestei Licente din diversele documente cu o singura copie care sa fie inclusa în colectie cu conditia sa urmati regulile acestei Licente pentru copii identice pentru fiecare document în toate celelalte privinte Puteti sa extrageti un document dintr-o astfel de colectie si sa-l distribuiti individual sub aceasta Licenta cu conditia de a include o copie a acestei Licente în documentul extras si sa urmati conditiile acestei Licente în toate celelalte privinte în legatura cu copiile identice ale acelui document 7 AGREGAREA CU LUCRARI INDEPENDENTE O compilatie a Documentului sau a unui derivat al sau cu orice document sau lucrare separata independenta, în sau pe un volum de stocare sau distribuire se numeste "agregat" daca drepturile de autor rezultate în urma compilarii nu sunt folosite pentru a limita drepturile legale ale utilizatorilor compilatiei mai mult decât permit lucrarile individuale Când Documentul este inclus într-un agregat, aceasta Licenta nu se aplica celorlalte lucrari din agregat care nu sunt ele însele rezultate derivate ale Documentului Daca cerintele legate de Textele de Coperta din sectiunea 3 se aplica acestor copii ale Documentului, atunci daca Documentul este mai putin de jumatate din întregul agregat atunci Textele de Coperta ale Documentului pot fi puse pe coperte care sa separe Documentul în cadrul agregatului, sau pe un echivalent electronic al acestora, daca Documentul se prezinta în format electronic Altfel ele trebuie sa apara pe copertele tiparite care îmbraca întreg agregatul 8 TRADUCERE Traducerea este considerata o forma de modificare, drept care puteti distribui traduceri ale Documentului sub cerintele sectiunii 4 Înlocuirea Sectiunilor Invariante cu traduceri ale acestora necesita permisiune speciala din partea celor care detin drepturile de autor, însa puteti include traduceri ale unora dintre sau tuturor Sectiunilor Invariante împreuna cu variantele originale ale acestora Puteti include o traducere a acestei Licente cât si toate notificarile de licenta din Document, cât si Limitarile de Responsabilitate atâta timp cât includeti si versiunea originala în engleza a acestei Licente, plus versiunile originale ale respectivelor notificari de licenta si limitari de responsabilitate În cazul aparitiei oricaror discrepante între versiunea tradusa si versiunea originala a acestei Licente, a vreunei notificari de licenta sau a vreunei limitari de responsabilitate, versiunea originala are prioritate Daca vreo sectiune din Document este Numita "Multumiri", "Dedicatii" sau "Istorie" cerinta (din sectiunea 4) de a-i Pastra Titlul (sectiunea 1) va necesita în mod normal schimbarea titlului în sine 9 REZILIERE Nu puteti copia, modifica, sublicentia sau distribui Documentul decât în conditiile specificate explicit în aceasta Licenta Orice copiere, modificare sau redistribuire a Documentului în vreo alta conditie este nula si va va anula în mod automat drepturile conferite de aceasta Licenta Pe de alta parte, tertilor carora le veti fi transmis copii sau drepturi în conformitate cu aceasta Licenta nu li se vor anula aceste drepturi atâta timp cât i se conformeaza 10 VERSIUNI VIITOARE ALE ACESTEI LICENTE Fundatia Free Software (Free Software Foundation) poate publica versiuni noi, revizuite ale acestei Licente (GNU Free Documentation Licence) din timp în timp Aceste noi versiuni vor pastra spiritul acestei versiuni dar pot diferi în privinta detaliilor, cu scopul de a se adresa unor noi probleme reale sau potentiale Vezi http://www gnu org/copyleft/ Fiecarei versiuni ale acestei Licente îi este asociat un numar de versiune distinct Daca Documentul specifica un anumit numar de versiune "sau orice versiune ulterioara" al acestei Licente, aveti de ales între a va conforma termenilor si conditiilor ori ale versiunii specificate explicit sau ale oricarei variante ulterioare publicate (nu ca varianta preliminara) de catre Free Software Foundation Daca Documentul nu specifica un numar de versiune al acestei Licente atunci puteti alege orice versiune publicata (nu ca varianta preliminara) de catre Free Software Foundation INTRODUCERE - 7 - ÎNCOVOIEREA STATICA ÎNCOVOIEREA STATICA - 68 - - 69 - ECUATII, ALGORITMI SI SOLUTII ECUATII, ALGORITMI SI SOLUTII EVALUARI NUMERICE EVALUARI NUMERICE VIBRATII DE ÎNCOVOIERE VIBRATII DE ÎNCOVOIERE VIBRATII INDUSE DE CURGEREA PARALELA VIBRATII INDUSE DE CURGEREA PARALELA CONCLUZII 